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In recent years, improved astrometric and spectroscopic techniques have opened the possibility of 
measuring the temporal change of radial and transverse position of sources in the sky over relatively 
short time intervals. This has made at least conceivable to establish a novel research domain, which 
we dub "real-time cosmology". We review for the first time most of the work already done in this 
field, analysing the theoretical framework as well as some foreseeable observational strategies and 
their capability to constrain models. We first focus on real-time measurements of the overall redshift 
drift and angular separation shift in distant sources, which allows the observer to trace the back- 
ground cosmic expansion and large scale anisotropy, respectively. We then examine the possibility 
of employing the same kind of observations to probe peculiar and proper accelerations in clustered 
systems, and therefore their gravitational potential. The last two sections are devoted to the future 
change of the cosmic microwave background on "short" time scales, as well as to the temporal shift 
of the temperature anisotropy power spectrum and maps. We conclude revisiting in this context 
the usefulness of upcoming experiments (like CODEX and Gaia) for real-time observations. 
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I. INTRODUCTION 

In 1920, Willem de Sitter made this remark: "The choice between the systems A and B is purely a matter of taste. 
There is no physical criterion as yet available to decide between them". He referred to system A and B, respectively, as 
the Einstein's solution to a quasi-static Universe filled with both matter and a cosmological constant A (the latter not 
yet considered as a source term in the Einstein equations, but just as an additional geometrical term which preserved 
covariance) and his solution to an "empty" universe with A. Almost 100 years after de Sitter's comment, cosmology 
has entered the so-called "era of precision": the content of the Universe has been established - at least with respect to 
the gravitational properties of matter - with very high confidence. 

The overall scenario where we, as observers, are supposed to live in turned out to be a universe with nearly flat 
geometry on large scales, composed by a bit less than 30% of matter and a bit more than 70% of a nearly uniform 
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component (called dark energy) with pressure negative enough to drive a late time accelerated expansion. Only a 
small percentage of matter (~ 4%) is made of baryons, while the remaining amount (~ 20%, called dark matter) does 
not directly emit electromagnetic radiation. The leading idea behind the formation of cosmic structures is that they 
arose by gravitational instability from small perturbations in the distribution of dark matter in the early Universe, 
eventually dragging baryons in the potential wells. This scenario is often called "concordance model", since it is 
consistent with the vast majority of cosmological observations to date. 

The leading cosmological observable is, at present, provided by measurements of the Cosmic Microwave Background 
(CMB) anistropies: small perturbations in the baryon-photon plasma in the early Universe left an imprint on the 
photon temperature and polarization anisotropy pattern at the time of decoupling, resulting in a goldmine of crucial 
information for cosmologists. In particular, the radiation we see as the CMB appears to come from a spherical surface 
around the observer such that the radius of the shell is the distance each photon has traveled since it was last scattered 
at around the epoch of decoupling. The physical scale of the first acoustic peak of the CMB power spectrum is set by 
the sound horizon at recombination, while its angular scale uncovers the angular diameter distance to last scattering 
surface, at z ~ 1100 (see [1] for an introductory review on the subject). 

No less relevant are the recent measurements of baryon acoustic oscillations (BAO). The same scales related to the 
CMB anisotropy can be spot also in the baryonic matter distribution, leading to acoustic oscillations in the matter 
power spectrum. The main benefit of BAO comes from the concept of standard ruler: we can infer the distance 
of an object of known size measuring its angular dimension. In this case the known size is precisely predicted by 
linear theory of perturbations and corresponds to the size of the sound horizon at decoupling, measured by the CMB. 
Measuring the radial and tangential size of BAO at a certain redshift one is able to reconstruct the Hubble function 
H(z) and the angular diameter distance at that redshift (consult Ref. [2] for a detailed review on BAO). 

Last but not least, cosmologists recently learned to exploit the measurement of luminosity distance of type la 
Supernovae (SNIa) as a function of redshift. After the discovery of a specific correlation between the shape and the 
height of their light curve, SNIa started to be used as standard candles. The further SNIa appeared dimmer then one 
would have expected in a Einstein-de-Sitter model, that is in a Universe filled only with matter. This experiments 
unveiled the presence of a new form of energy (dark energy) that does not emit electromagnetic radiation and whose 
pressure is so negative that speeds up the expansion of the universe whenever it dominates the density budget (a 
recent review on the subject can be found in [3]). 

All these complementary and pivotal measurements (in addition to determination of the Hubble constant [4]) 
combined together to shape the previously mentioned concordance scenario. Also, weak lensing observations have 
recently started to play a prominent role (see [5]). However, these observables offer an information that is primarily 
focused on distances (angular diameter distance and luminosity distance) and perturbation growth. Reconstructing 
in an accurate way the background expansion means reconstructing in an accurate way H(z). This encodes a great 
part of the story of cosmological parameters (like the energy density and equation of state parameters, ili(z) and 
Wi(z) respectively) and hence the kinematics and dynamics of expansion. Unfortunately, distances are integrals over 
H(z), which in turn is an integral over Wi(z): therefore, the error propagation of these parameters is quite complex. 
In other words, the correspondence between distances and Wi(z) is not straightforward. In addition, degeneracies in 
parameter space and systematic errors might not guarantee that accuracies on the estimates significantly improve 
with time. Hence the advent of further observables able to test the expansion in a different way and in complementary 
redshift windows is becoming more and more appealing. 

In this context, given the advancement in technology occurred over the last forty years, the idea of measuring 
temporal variation of astrophysical observable quantities over a few decades, i.e. in real time, has become conceivable: 
hence, the possibility of establishing a novel and appealing research field, which we dubbed "real-time cosmology [6]. 
Real-time observations may be related to variations of radial and transverse position and/or velocity of a given source. 
In this paper we give, for the first time, an overview of most of the work already done on the topic. In particular, 
we deemed useful to divide real-time observables into two classes: temporal shifts mainly tracing the background 
expansion and temporal shift caused by peculiar motions. For each class one can think about two ideal sub-classes, 
corresponding to drifts in radial and transverse directions with respect to the observer line of sight. 

We will start by introducing the first and more ancient real time observable, namely the redshift drift. First 
conjectured by Sandage in 1962 [7], the redshift drift belongs to the first class of real-time observations and regards 
the temporal variation of redshift of distant sources as a tracer of the background cosmological expansion (i.e. in 
the radial direction). As we will show in this review, it could be a very important cosmological probe, since it is a 
straightforward measure of the change of H (z) with respect to its present value, and would provide a direct detection 
of acceleration in the expansion. Today, spectroscopy has reached already a sensitivity of a few meters per second. 
Lyman-a clouds along the line of sight of very distant and bright sources like quasars [8] are poorly affected by peculiar 
motions. By using spectra populated by many sharp lines, some authors have shown that a statistical sensitivity of 
the required level might already be reached with the next generation of optical telescopes [9] . 

By using the same real-time observable and selecting sources for which the cosmological background signal is 
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expected to be tiny (e.g. sources closer to us) it would also be possible to track the variation of peculiar velocity 
of objects in clusters and galaxies over a few decades, allowing the measurement of the acceleration caused by 
potential wells (which we call peculiar acceleration [10]). This procedure opens up the possibility of reconstructing 
the gravitational potential in a direct way, using this detection to distinguish between different gravity models, like 
for example Newtonian dynamics and the MOND paradigm [11]. 

Also belonging to the first class of real-time observables is the analogue of redshift drift, but in the transverse 
direction. In particular it has been shown that the temporal change of the angular separation between distant sources 
(e.g. quasars) can be used to detect a background anisotropic expansion. This is the so-called " cosmic paralla^ [6]. 
The standard model of cosmology rests on two main assumptions: general relativity and a homogeneous and isotropic 
metric, the Friedmann-Robertson- Walker metric (henceforth FRW). While general relativity has been tested with 
great precision at least in laboratory and in the solar system, the issue of large-scale deviations from homogeneity 
and isotropy is much less settled. There is by now abundant literature on possible tests of the FRW metric, and 
on alternative models invoked to explain the accelerated expansion by the effect of strong, large-scale deviations 
from homogeneity (see [12] for a review). In a FRW isotropic expansion the angular separation between sources is 
constant in time (except for the effect of peculiar motions) and the cosmic parallax vanishes. While deviations from 
homogeneity and isotropy are constrained to be very small from cosmological observations, these usually assume the 
non-existence of anisotropic sources in the late universe. Conversely, dark energy with anisotropic pressure may act as 
a late-time source of anisotropy. Even if one considers no anisotropic pressure fields, small departures from isotropy 
cannot be excluded, and it is interesting to devise possible strategies to detect them. The anisotropic expansion 
can be either intrinsically set by the metric itself (like in Bianchi models) or mimicked by an off-center position of 
the observer in a inhomogeneous Universe (like for example in Lemaitre-Tolman-Bondi (LTB) void models). In the 
latter case a detection of cosmic parallax would also provide a test of the Copernican Principle stating that we, as 
observers, do not occupy a favourite position in the Universe. If again the signal is dominated by peculiar motions in 
a bound system, like in our own Galaxy, the angular temporal shift becomes proportional to peculiar acceleration in 
the transverse direction, which we dub proper acceleration. 

This review is organized as follows. In the first two sections we focus on the first class of real-time observables, i.e. 
the ones mapping out the background expansion. Speficifically, in Section II we introduce the concept of redshift drift, 
its derivation and forecasted constraining power for several dark energy FRW models, as well as for less symmetric 
space-times. Then, in Section III, we examine the cosmic parallax signal in LTB models with an off-center observer, 
both from an analytical point of view and with a full numerical derivation. We present the cosmic parallax in Bianchi I 
models as well, where the source of anisotropy is an anisotropically distributed dark energy density. In the subsequent 
sections we concentrate on works that examine the second class of real-time observables, that is on signals generated by 
peculiar motions. In Section IV we derive the real-time peculiar acceleration expression both in linear approximation 
and in non-linear structures like clusters and galaxies. In Section V proper acceleration detected by temporal drift 
in the angular position of test particles in our own Galaxy is presented for the first time. Section VI is dedicated to 
recent papers presenting forecasts of future time variations of CMB temperature, angular power spectrum and maps 
[13, 14]. Finally, some details about the observational strategies, instrumental required accuracies (both astrometric 
and spectroscopic) and capability of already planned mission to measure real-time observables are collected in Section 
VII. In Section VIII we draw our conclusions. 



II. THE REDSHIFT DRIFT 

The measurement of the expansion rate of the Universe at different redshifts is crucial to investigate the cause of 
the accelerated expansion, and to discriminate candidate models. Until now, a number of cosmological tools have 
been successfully used to probe the expansion and the geometry of the Universe. 

Depending on the underlying cosmological model, one expects the redshift of any given object to exhibit a specific 
variation in time, or redshift drift. An interesting issue, then, is to study whether the observation of this variation, 
performed over a given time interval, could provide useful information on the physical mechanism responsible for the 
acceleration, and be able to constrain specific models. This is the main goal of this Section. In addition to being a 
direct probe of the dynamics of the expansion, the method has the advantage of not relying on a determination of the 
absolute luminosity of the observed source, but only on the identification of stable spectral lines, thus reducing the 
uncertainties from systematic or evolutionary effects. 

The possibility of using the time variation of the redshift of a source as probe of cosmological models was first 
proposed by Sandage [7]. The predicted signal was less than a cm/s per year and, at the time, deemed impossible 
to observe. Although the test was mentioned a few other times in the literature over the last decades (e.g. [15-17]) 
it was not until recently that the feasibility of its observation was reassessed by Loeb [8] and judged within the 
scope of future technology (see also [18]). In particular, the foreseen development of extremely large observatories, 
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such as the European ELT (E-ELT), the Thirty Meter Telescope (TMT) and the Giant Magellan Telescope (GMT), 
with diameters in the range 25-100 m, and the availability of ultra-stable, high-resolution spectrographs, encouraged 
new evaluations of the expected signal for the current standard cosmological model (dominated by a cosmological 
constant), through the analysis of realistic simulations. The conclusion of such studies was that the perspective for the 
future observation of redshift variations looks quite promising. For example, the authors in [19] pointed out that the 
CODEX (COsmic Dynamics Experiment) spectrograph should have the right accuracy to detect the expected signal 
by monitoring the shift of Lyman-a forest absorption lines of distant (z > 2) quasars over a period of a few decades. 
These sources have the advantage of being very stable and basically immune from peculiar motions. In Section VII A 
we will explore in more detail the observational strategy. 

Such new prospects lately prompted renewed interest in the theoretical predictions of the redshift variation in 
different scenarios (see e.g. [20-29]). 

Despite its inherent difficulties, the method has many interesting advantages. One is that it is a direct probe of the 
dynamics of the expansion, while other tools (e.g. those based on the luminosity distance) are essentially geometrical 
in nature. This could shed some light on the physical mechanism driving the acceleration. For example, even if the 
accuracy of future measurements will turn out to be insufficient to discriminate among specific models, this test would 
be still valuable as a tool to support the accelerated expansion in an independent way, or to check the dynamical 
behaviour of the expansion expected in general relativity compared to alternative scenarios. It must be noted that 
radial BAO surveys can also be used to measure H(z). This is due to the fact that radial BAO are a measure of 
the comoving distance in a given redshift bin, which, for a narrow enough bin, is inversely proportional to -ff (zbin)- 
Nevertheless, BAO measurements are non-trivial and are subject to their own systematics. On the other hand the 
redshift drift, despite being observationally challenging, is conceptually extremely simple. For example, it does not 
rely on the calibration of standard candles (as it is the case of type la SNe) or on a standard ruler which originates 
from the growth of perturbations (such as the acoustic scale for the CMB) or on effects that depend on the clustering 
of matter (except on scales where peculiar accelerations start to play a significant role). Therefore, the redshift drift 
will also serve as a useful cross-check for radial BAO. Third, by using distant quasars, it will provide constraints on 
the cosmic expansion at redshifts z > 2, where supernovae and large scale surveys have difficulties in providing quality 
data. Finally, it allows to distinguish between true acceleration, as for dark energy models, and apparent acceleration, 
as in void models, as we will discuss in Section II B. 



A. The redshift drift in homogeneous and isotropic universes 

The basic theory behind redshift variation in time is quite simple. One starts assuming that the metric of the 
Universe is described by the simplest FRW metric. The observed redshift of a given source, which emitted its light 
at a time t s , is, today (i.e. at time to), 

*.(*>) = ^-i, a) 

a(t a ) 

and it becomes, after a time interval Ato (At s for the source) 

z a (t + At ) = a ^±^\-l. (2) 
a{t s + At s ) 

The observed redshift variation of the source is, then, 

^ _ a(t + At ) _ a (t ) 
a(t s + At s ) aits)' 

which can be re-expressed, after an expansion at first order in At/t, as: 



(3) 



Clearly, the observable Az s is directly related to a change in the expansion rate during the evolution of the Universe, 
i.e. to its acceleration or deceleration, and it is then a direct probe of the dynamics of the expansion. It vanishes if 
the Universe is coasting during a given time interval (i.e. neither accelerating nor decelerating). We can rewrite the 
last expression in terms of the Hubble parameter H(z) — a(z)/a(z): 

Az s = H At (l + zs- , (5) 
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where we have dropped the subscript for simplicity. The function H(z) contains all the details of the cosmological 
model under investigation. Finally, the redshift variation can also be expressed in terms of an apparent velocity shift 
of the source, Av = cAz s /(l + z s ). 
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Table I. Expansion rate for several dark energy models in the framework of homogeneous and isotropic cosmologies. The 
redshift drift evolution corresponding to these Hubble functions is depicted in Fig. 1. In the affine model Cl m = (po — Pa)/ Peru 
(for more detailed designation of all the parameters we refer to [22]). 

In [22] the authors analysed many currently viable dark energy cosmological models based on the assumption of 
homogeneity and isotropy. Their Hubble expansion rates as a function of the scale factor are collected in Table I. 
These models have often been invoked as candidates to explain the observed acceleration [30], i.e. they have not been 
falsified by available tests of the background cosmology. For each class the best fit values found in [30] was assumed, 
and the parameters were varied within their 2a uncertainties. Clearly some models may be preferred with respect to 
others based on the fact that they fit the data well with a smaller number of parameters. Nevertheless, it is interesting 
to explore as many models as possible, since future observations of the time variation of redshift could reach a level of 
accuracy which could allow to better discriminate competing candidates, and to understand the physical mechanism 
driving the expansion. 

In order to perform a forecast analysis the predicted accuracy of observations expected from an experiment like 
CODEX (see Section VII A) was adopted. The latter is entirely based on the Monte Carlo simulations and discussed 
by [19]. As expected, in its simplest formulation, the accuracy scales as the square root of the total number of quasars 
and is a decreasing function of redshift: the more the farther the source is from us observer the higher is the number 
of features captured in the Lyman-a spectra. Details about the experimental accuracy and discussions related to it 
can be found in Section VII A. 

These error bars were used to construct simulated data and get a feeling of the possible constraints to viable dark 
energy models. Fig 1 shows the predicted signal for the models assembled in Table I. All the predictions were derived 
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assuming At = 30 years and a future dataset containing a total of 40 quasars spectra uniformly distributed over 5 
equally spaced redshift bins in the redshift range 2-5 with a S/N=3000, observed twice over the aforementioned time 
span. This observational strategy was properly justified in [22] (details on this assumptions will be encountered again 
in Sec. VII A). 

From Eq. 5 it is clear that the expected velocity shift signal increases linearly with At, so that it is straightforward 
to calculate the expected signal when a different period of observation is assumed. It is clear that the observation 
of velocity shift alone can be affected by the degeneracies of the parameters that enter H(z), limiting its ability to 
constrain cosmological models. The uncertainties on parameter reconstruction (particularly for non-standard dark 
energy models with many parameters) can be rather large unless strong external priors are assumed. When combined 
with external inputs, however, the time evolution of redshift could discriminate among otherwise indistinguishable 
models. 

In Ref. [22] a Fisher matrix analysis allowed to estimate the best possible accuracy attainable on the determination 
of the parameters of a certain model. Given a set of cosmological parameters p,-, i = 1, ...,n, and the corresponding 
Fisher matrix Fij (that is easily calculated based on a theoretical fiducial model and the assumed data errors), the 

1 11 

best possible la error on pi is given by Api = CJ , where the covariance matrix CV, is simply the inverse of the 
Fisher matrix: Cy = F^j . The prospect of detecting departures from the standard ACDM case could in principle be 
one of the real assets of observing the time evolution of redshift, and is thus worthy of closer investigation. Since the 
simulated data used in the analysis assume that quasars are used as a tracer of the redshift evolution, we expect that 
the more constrained models will be those that have the largest variability in the redshift range 2 < z < 5. It is at 
least conceivable that suitable sources at lower redshifts than those considered in this work could be used to monitor 
the velocity shift in the future. This would be extremely valuable, since some non-standard models have a stronger 
parameter dependence at low and intermediate redshifts (see Fig. 1), that could be exploited as a discriminating tool. 
In [8] speculative possibilities of using other sources have been indicated, like masers in galactic nuclei, extragalactic 
pulsars or gravitationally lensed galaxy surveys: this would further extend the lever arm in redshift space and increase 
the ability of constraining models. These may certainly be interesting topics for further studies. 

Assuming that the fiducial model has tt\ — 0.7 and = and that both J7a and fi/c can vary, Ref. [22] found 
A£1a = 0.2 and ASlj. = 0.25 at la. Fixing = 0, the bound on fl\ becomes 0.007 (at ltr). If dark energy is modelled 
by a constant equation of state (with a fiducial value w = — 1) and the flatness constraint is imposed we find a looser 
bound on the dark energy density, Aflde — 0.016, and quite a large error on the equation of state, Aw = 0.58. This 
clearly shows that different assumptions on the knowledge of any parameter has an influence on all the others. The 
parameter fi^ can be much better constrained using external datasets, such as the CMB anisotropy. 

The DGP model is the one for which the tightest constraints have been obtained: A£l r = 0.0027 at la, assuming 
il r = 0.13 as a fiducial value. This is not only due to the strong dependence of the velocity shift on fi r (see Fig. 1), 
but also to the simplicity of the model, which depends on only one parameter (in this respect, this is the simplest 
model, together with the standard flat ACDM). In general, it is to be expected that models with less parameters 
perform better. 

For what concerns the Chaplygin model, when both A and 7 vary freely, no interesting constraint can be obtained 
observing the velocity shift with the assumed QSO data: we find AA = 0.42 and A7 =1.4 (for the fiducial values 
A = 0.7 and 7 = 0.2). Fixing A, on the contrary, results in a very tight bound on 7: A7 = 0.008. 

The interacting dark energy and the affine equation of state models show a large variability in the redshift range 
we are exploring: this is to be expected, since in both models the matter-like component departs from the usual 
a~ 3 scaling, giving a distinctive signature when one looks at higher and higher redshifts. If is known, the affine 
parameter a can be reconstructed with an error Aa = 0.005. If, in addition, w = —1 we find A£ = 0.06 for the 
interacting dark energy model (for a fiducial value £ = 3). 

The other models do not seem to have very interesting signatures to be exploited, at least in the redshift range 
considered in our analysis. 

Fig. 2 shows the comparison among the predicted velocity shifts for all the models described earlier, assuming 
parameter values that are a good fit to current cosmological observations (including the peak position of the CMB 
anisotropy spectrum, the SNe la luminosity distance, and the baryon acoustic oscillations in the matter power spec- 
trum). In other words, the models shown in Fig. 2 cannot be easily discriminated using current cosmological tests of 
the background expansion. If we assume that the ACDM model is the correct one, and simulate the corresponding 
data points for the velocity shift, using a \ 2 test we can quantify how well we can exclude the competing models 
based on their expected signal. As it is clear from Fig. 2, some models can be excluded with a high confidence level. 
In particular, the Chaplygin gas model and the interacting dark energy model would be excluded at more than 99% 
confidence level, and that the affine model would be out of the 1 a region. 

All the results above were obtained assuming an equal number (8) of quasars for each of 5 redshift bins centered at 
z = {2, 2.75, 3.5, 4.25, 5}, all of the same redshift width of 0.75 (such a uniform distribution was also assumed in the 
simulations performed by [19]). A good estimate of the error bars for such binning can be achieved using the brightest 
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quasars from the SDSS DR7 as a benchmark [25]. This is explained in more details in Appendix B. Moreover, one 
should notice that the effective S/N on the measured data points decrease as At 3 / 2 , as the signal is linear in time and 
the noise scales as Ai -1 / 2 , and can become significantly higher over only a few decades of observations. 



B. The redshift drift in less symmetric space-times 



FRW universes are based on a maximally symmetric space-time derived from the assumption of the cosmological 
principle: the observable universe is homogeneous and isotropic (in the weaker and matter-of-fact version at least on 
large scales). The cosmological principle itself is originated by a combination of the Corpernican principle, stating 
that we are not located at a favoured position in space, and the observed isotropy [31]. The inferred existence of a 
late-time accelerated expansion (and consequently the hypothesis of a dark energy component) has its main foundation 
on the conjectured flatness of space-time justified by CMB observations together with the Copernican principle. The 
verification of the latter has recently attracted attention to test whether the relaxation of its statement may help 
explain recent time expansion without invoking the existence of a new exotic component or modified gravity theory. 

In this framework, less symmetric space-times have been considered; in particular, since one should expect to observe 
the redshift drift in any expanding space-time, in [32] the authors have derived an expression for this observable in 
spherically symmetric universes, where the observer is located at the centre. 

One type of possibly useful coordinates that allow for a general treatment are the observable coordinates {w, y, 9, ip}, 
where w marks the past light-cones of events along the worldline C of the observer. The metric is 

ds 2 = -A 2 (w, y)dw 2 + 2A(w, y)B(w, y)dydw + C 2 {w, y)dSl 2 , (6) 

spherically symmetric around the world-line C defined by y = 0. The redshift is given by 

^ _|_ ^ (u a ka) emission A(wq,0) 

(u Q ,fc a ) bserver A{Wo,y) 

for a given value Wq. Here the matter velocity and photon wave-vector are u a = A~ 1 5'^ and k a = (AB)~ 1 5y , 
respectively. The isotropic expansion rate around an observer located at the centre is defined as H = V a u"/3 and 
eventually turns out to be: 

rj( s 1 ( d w B{w,y) d w C(w,y) ^ 

In the case of a dust-dominated universe the covariant derivative of u a takes the form V a up = H(g a p + u a up) + a a p, 
where a a p is the traceless and symmetric shear satisfying u Q er Q( 3 = 0. 
The expression for the redshift drift is then derived from Eq. (7) 

Sz(w ,y) = ( d w A(w 0) 0) d w A(w ,y) \ 

Sw 1 >\ A{w o ,0) ~ A(w ,y) )' [ ' 

Choosing A(wq,0) = 1 and y such that d w B(wo,y) — d w A(wo,y), it follows that 

5z(w , y) = + _ _ 1 

5w V3 

where a — a a ^a a ^j2 is the scalar shear. Indeed this is the general form for the redshift drift measured by an observer 
located at the centre of a spherically symmetric universe. In the following, we will present numerical and analytical 
calculations of the redshift drift in some LTB models also including off-centre observers. 



1. The redshift drift in Lemaitre-Tolman-Bondi void models 



The most general metric describing a spherically symmetric and inhomogeneous universe in comoving coordinates 
is the LTB metric. In particular, the latter is well designed to describe universes where the observer is located near 
the centre of a large void embedded in an Einstein-de Sitter cosmology (see e.g. [33-36]). Models with void sizes 
which, although huge by any means, are "small" enough (z ~ 0.3 — 0.4) not to be ruled out due to distortions of the 
CMB blackbody radiation spectrum [37] are capable of fitting the observed SNIa Hubble diagram and the CMB first 
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Figure 1. The apparent spectroscopic velocity shift over a period Ato = 30 years, for a source at redshift z 3 , for the models 
described in Table I. From Ref. [22]. 
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Figure 2. The predicted velocity shift for the models presented in Table I, compared to simulated data as expected from the 
CODEX experiment. The simulated data points and error bars are estimated from Eq. (103), assuming as a fiducial model the 
standard ACDM model. The other curves are obtained assuming, for each non-standard dark energy model, the parameters 
which best fit current cosmological observations. From Ref. [22]. 



peak position and compatible with the COBE results of the CMB dipole anisotropy, as long as the observer is not 
too far from the center [25, 38] (although recent works [39-41] claim that the blackbody spectrum of the CMB, in 
combination with the kinetic Sunyaev-Zel'dovich effect due to free electrons effectively rule out these models). The 
off-center displacement is limited to ~ 200 Mpc by supernovae, whereas the CMB dipole limit it to somewhere in 
between 30 and 60 Mpc depending on some a priori assumptions [25, 38, 42]. 

In principle, the redshift drift in such models needs an exact treatment where the full relativistic propagation of 
light rays is taken into account. We will begin by introducing the Einstein equations in such a metric and we will 
then present the light geodesic equations. 

The LTB metric can be written as (primes and dots refer to partial space and time derivatives, respectively): 



ds 2 



-dt 2 



[R'(t,', 



■dr 2 + R 2 {t 1 r)dQ 2 1 



(11) 



where j3{r) can be loosely thought of as a position dependent spatial curvature term. Two distinct Hubble parameters 
corresponding to the radial and perpendicular directions of expansion are defined as 



H\\ = 



R'/R' , 
= R/R. 



(12) 
(13) 



Note that in a FRW metric R — ra(t) and H\\ = H±. This class of models exhibits implicit analytic solutions of the 
Einstein equations in the case of a matter-dominated universe, to wit (in terms of a parameter rf) 



R = (cosh?? - 1)— + i?iss 
V/3t = (sinhry - r?) ^ + i?i s; 



, a + /3R lss . 
cosh rj + \l — — smli 7] 

pR\ss 



sinh rj 



la + PR 



Iss 



PR 



Iss 



(cOSh 7] — 1) 



(14) 
(15) 



where a, j3 and R\ ss are all functions of r. In fact, R\ ss (r) stands for R(0,r) and we will choose t = to correspond 
to the time of last scattering, while a(r) is an arbitrary function and (3(r) is assumed to be positive. In performing 
the calculations, it is much simpler to set Ri ss = 0, as i?i ss ^ introduces arduous numerical problems: one only has 
to note that they will not be valid for z > 10, as has to be the case anyway. 
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Figure 3. The annual redshift drift for different models assuming an observer at the center. The upper, blue solid lines represent 
the ACDM model. The green, dashed line corresponds to a self-accelerating DGP model with Q. Tc — 0.13. The dot-dashed 
lines stand for the 3 void models considered here: the dark brown (indistinguishable) lines are for Models I and II, while the 
red line just above correspond to the cGBH model. The line at the bottom of the plot corresponds to an universe with only 
matter in a FRW metric (the CDM model). From Ref. [25]. 



The full geodesic equations are written in Appendix A, where an algorithm is provided to compute both the redshift 
drift and the cosmic parallax (see next Section). The former effect was first calculated in [6] (see Sec. Ill)) for two 
distinct specific LTB models, while the latter was thoroughly investigated in [25] for the same models as well as for a 
third one. 

Two of these three models (Ref. [38, 43]) are characterized by a smooth transition between an inner void and an 
outer region with higher matter density and described by the functions: 



a(r) = (Hr*) 



Aa I , r — r vn 

1 1-tanh 



2 V 2Ar 



(16) 



Kr) = (Hl%)\ 2 ^ (l - tanh r -^j , (17) 

where Aa, r vo and Ar are three free parameters and -ff° u g is the Hubble constant at the outer region, set at 
51 km s _1 Mpc -1 . We will dub the two models "Model I" and "Model II", and define them by the sets {Aa = 
0.9, r vo = 1.46 Gpc, Ar = 0.4r vo } and {Aa = 0.78, r vo = 1.83 Gpc, Ar = 0.03r vo }, respectively. These values of 
r vo correspond, in physical distances (that is, the distance obtained by multiplying the comoving separation by the 
scale factor), to void sizes of 1.34 and 1.68 Gpc, respectively. The third model under consideration is the so-called 
"constrained" model proposed in [44] referred to as the "cGBH" model. For this model, the parameters were chosen to 
maximize the likelihood as obtained in [44]; this model can be written in terms of a and /? ([25]). The main difference 
between the three models is that Model II features a much sharper transition from the void and that the cGBH model 
is almost twice as large. Although in principle the redshift drift will depend on the source position in the sky for 
off-center observers, it was shown in [25] that, unless the LTB models violate the CMB dipole measurements, the 
differences across the sky are less than 5%. 

The 3 LTB models we shall investigate in this and in the following section are depicted in Fig. 4. 

Fig. 3 illustrates the redshift drift as a function of redshift for ACDM, the DGP model, the old matter dominated 
model (CDM) and the 3 different void models. As could be expected, the void models predict a curve which is in 
between CDM and ACDM. 

As can be seen from Fig. 2, in many dark energy models the redshift drift is positive at small redshift, but becomes 
negative for z > 2. On the other hand, a giant void mimicking dark energy produces a very distinct z dependance 
of this drift, and in fact one has that dz/dt is always negative (see Fig. 3). This property and its potential as 
discriminator between LTB voids and ACDM was first pointed out in [45]. Moreover, it was recently proven in [29] 
that under certain general conditions the property dz/dt < always holds in void models. 

Fig. 5 depicts the redshift drift for both 7 and 13 years time span ([25]). Also plotted are the forecasted error 
bars obtainable by CODEX at the E-ELT. The correspondent possible achievable significance of detection is listed in 
Table II. As can be seen from this table, void models can be ruled out with a ~ 5<r confidence-level within a decade 
of observations. 
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Figure 4. Q m o and hii = H^/Ho for Model I (solid), Model II (dashed curve) and the cGBH model (red, long-dashed) as a 
function of the physical (as opposed to comoving) distance X. Note that the definition for Q m o we use differ from the one 
in [38, 43] and we do not get their characteristic over-density bump (or shell) surrounding the void. 




Figure 5. Redshift drift for different dark energy models for a total mission duration of 7 (left) and 13 (right) years and 
CODEX forecast error bars. In each plot, the upper 3, solid lines represent const. wCDM models for w — —1.25 (uppermost), 
w = — 1 (second) and w = —0.75 (third uppermost). The green, dashed line corresponds to a self-accelerating DGP model 
with Q,. c = 0.13. The three dot-dashed lines at the bottom of the plot stand for the 3 void models considered here: the dark 
brown (indistinguishable) lines are for Models I and II, while the red line just above correspond to the cGBH model. Similar 
to figures from Ref. [25], except for the inclusion of horizontal error bars. 



The predicted accuracy assumed in [25] although adopting the same expression as in Eq. (103) (except for the 
redshift power law exponent, which is not substantially relevant), leads to an error bar that is somewhat higher at 
higher redshift. This is the outcome of two factors: quasars which are at high redshift are usually dimmer, so their 
S/N is lower, and that may win over the power law factor. Moreover, the time spent observing each quasar is assumed 
to be the same, and this accounts for larger error bars at high redshift due to the lower apparent magnitudes of the 
corresponding quasars. On the other hand, in [22] the relative integration time for these sources was assumed to be 
increased in order to achieve the same average signal-to-noise ratio at all redshift bins. The observational strategy 
will be better analyzed in Sec. VII. In order to estimate the error bars in the left panel of Fig. 5, Ref. [25] made use 
of available SDSS quasars, selecting the brightest quasars in each redshift bin using the appropriate band for such 
bin. These details are covered in the Appendix B. In the right panel of the same figure a similar plot is shown from 
Ref. [28], where the authors also included redshift errors. 



III. COSMIC PARALLAX 



In every anisotropic expansion the angular separation between any two sources varies in time (except for particular 
sources aligned on symmetry axes), thereby inducing a cosmic parallax effect [6[. This is totally analogous to the 
classical stellar parallax, except here the parallax is induced by a differential cosmic expansion rather than by the 
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Model 


7 years 


10 years 


13 years 


Models I / II 


3.0a 


6.2a 


10a 


cGBH Model 


2.2a 


4.3a 


7.1a 



Table II. Estimated achievable confidence levels by the CODEX mission in 7, 10 and 13 years. 



observer's own movement. In other words, cosmic parallax will be present whenever one has shear in cosmology. An 
anisotropic expansion can either be experienced by an off-centre observer in an inhomogeneous and isotropic universe 
or a centered observer embedded in an intrinsically anisotropic expansion. LTB void models belong to the former 
class, and cosmic parallax was extensively analysed in [6, 25, 46], while Bianchi I models are specific example of the 
latter class and were studied in [47]. 

Measuring the cosmic parallax is therefore an independent test of late-time anisotropics in the Universe. More 
rigorously, it is a test of late-time anisotropic expansion, or shear. Competing tests include direct reconstruction of 
the Hubble diagrams using Supernovae [42, 43] and the CMB dipole [48] and quadrupole [49]. 



A. Cosmic parallax in Lemaitre-Tolman-Bondi void models 

As discussed in Section II B 1, LTB universes have two different Hubble parameters and therefore appear anisotropic 
to any observer except the central one. Therefore any such observer will see a sky affected by cosmic parallax. In 
fact, the amount of anisotropy is at first order directly proportional to this off-center distance. In void models, the 
cosmic parallax is also an independent test of the Copernican principle but without the degeneracy with our peculiar 
velocity that afflicts the constraints from the cosmic microwave background dipole [37, 38, 48]. This fact was used 
in [6, 25] to evaluate the feasibility of detecting the cosmic parallax in dark energy motivated void models. 

Due to their large distances and point-like properties, quasars are the obvious choice for observing the cosmic 
parallax. 

Figure 6 depicts the overall scheme describing a possible time-variation of the position of a pair of sources that 
expand radially with respect to the center but anisotropically with respect to the observer. We label the two sources a 
and &, and the two observation times 1 and 2. In what follows, we will refer to (t. r. 9, <j)) as the comoving coordinates 
with origin on the center of a spherically symmetric model. Peculiar velocities aside, the symmetry of such a model 
forces objects to expand radially outwards, keeping r, 9 and (f> constant. 



1. Estimating the parallax 

Following Fig. 6, let us assume an expansion in a flat FRW space from a "center" C observed by an off-center observer 
O at a distance X Q i, s from C . Since we are assuming FRW it is clear that any point in space could be considered a 
"center" of expansion: it is only when we will consider a LTB universe that the center acquires an absolute meaning. 
The relation between the observer line-of-sight angle £ and the coordinates of a source located at a physical radial 
distance X (corresponding to a comoving radial distance r) and angle 9 in the C-frame is 

p _ X cos 9 - A obs 

C ° S ^ (A2 + A 2 bs -2X obs Acosfl)V 2 ' [W) 

where all angles are measured with respect to the CO axis. We follow the approach of Ref. [6] and assume for 
simplicity (and clarity) that both sources share the same <fi coordinate. 

Consider first two sources at location oi, b± on the same plane that includes the CO axis with an angular separation 

71 as seen from O, both at distance X from C. After some time At, the sources move to positions 0.2,62 and the 
distances X and A b s will have increased by A t X and A t A b s respectively, so that the sources subtend an angle 

72 (see Fig. 6). In a FRW universe, these increments are such that they keep the overall separation 7 constant. 
However, if for a moment we allow ourselves the liberty of assigning to the scale factor a(t) and the H function a 
spatial dependence, a time-variation of 7 is induced. The variation 

A t 7 ee 71-72 (19) 

is the cosmic parallax effect and can be easily estimated if we suppose that the Hubble law is just generalized to 

A t X = XH(to, X)At ee XHxAt , (20) 



13 




Figure 6. Cosmic parallax in LTB models. C stands for the center of symmetry, O for the off-centre observer, a and b for 
two distinct distant light sources, such as quasars. For these latter three, the subindex 1 and 2 refer to two different times 
of observation. For clarity purposes we assumed here that the points C, 0,ai,bi all lie on the same plane. By symmetry, 
points a2,b2 remain on this plane as well. Comoving coordinates r and ro correspond to physical coordinates X and Xq. The 
difference between the angular separation of sources, At"/ = 71 — 72, is the cosmic parallax. The angular separation 74, in turn, 
is calculated as the difference between the angle £ of the incoming geodesies coming from a and b at time t (71 = £ a i — £6i). 
From Ref. [6]. 



where 



X(r) = J gU 2 dr' = J a(t ,r')dr' , (21) 

generalizes the FRW relation Afrw = a(^o) r i n a metric whose radial coefficient is g rr . 

For two sources a and b at distances much larger than X b s (which in practice in usual models corresponds to 
z a ,b ^ 0.1), after straightforward geometry we arrive at 



A t7 = AtX ohs 



{H hs — H a ) — (.Hobs — Hb) : ' 

X a Xb 



(22) 



It is important to note that this simple analytical estimate have been verified numerically, and the angular dependence 
of the cosmic parallax for sources at similar distances has been verified to hold to very high precision. As can be seen 
above, the signal A t 7 in (22) depends both on the sources' positions on the sky (the angles 9 a ,b) an d on their radial 
distances to the center (X a and Xb). In what follows we will consider two simplified scenarios for which the sources 
lie either: (i) on approximately the same redshift but different positions; (ii) on approximately the same line-of-sight 
but different redshifts. 

For case (i), we can average over 9 a .b to obtain the average cosmic parallax for two arbitrary sources in the sky 
(still assuming they lie on the same plane that contains CO). If both sources are at the similar redshifts z a ~ Zb = z 
(corresponding to a physical distance X), then the average cosmic parallax effect is given by 

(A t7 ) pcrp ~ ^| ^ |sin0 o -sin0 6 |d0 a d0b = At (H ohs - H x ) . (23) 

where we defined a convenient dimensionless parameter s such that 



»-obs 



x 



< 1 . (24) 



Note that at this order the difference between the observed angle £ and can be neglected [6]. We can also convert 
the above intervals AX into the redshift interval Az by using the relation r = J dz/H(z). Using (21) we can write 
AX = a(to, X)Az/H(z) <~ Az/H(z) (we impose the normalization a (to, A Q b s ) = 1), where H(z) = H(t(z),X). One 
should note that in a non-FRW metric, one has s 7^ ro/r. 

In a FRW metric, H does not depend on r and the parallax vanishes. On the other hand, any deviation from FRW 
entails such spatial dependence and the emergence of cosmic parallax, except possibly for special observers (such as 
the center of LTB). A constraint on A t 7 is therefore a constraint on cosmic anisotropy. 
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Rigorously, the use of the above equations is inconsistent outside a flat FRW scenario; one actually needs to perform 
a full integration of light-ray geodesies in the new metric. Nevertheless, following [6] we shall assume that for an order 
of magnitude estimate we can simply replace H with its space-dependent counterpart given by LTB models. In order 
for an alternative LTB cosmology to have any substantial effect (e.g., explaining the SNIa Hubble diagram) it is 
reasonable to assume a difference between the local Hobs and the distant Hx of order i? bs [38]. More precisely, 
putting H b s — Hx — H b s Ah then using (23) one has that the average A t "f is of order 



(A*7> 



20 s Ah /xas/year 



(25) 



for two sources at the same redshift. 

Similarly, for case (ii) we haver source pairs at approximately same position 8 but different (yet similar) redshifts, 
and one has (using (22)) 



A t7 



['ad 



Az 

s sinOAhAt Az/X /ias/year ~ 20 s sin 6 Ah — /xas/year, 



(26) 



where it was assumed that X ~ zH(z) l . The average radial cosmic parallax for sources between 10 and 200 times 
Xobs can be obtained numerically to be 



<A t 7> 



rad 



At (H ohs - H x )sm9 


/■200 /-200 


1 1 


190 2 


J 10 J 10 


Sa S b 



Therefore, one can estimate for the radial signal 

(A t7 ) 



rad 



d(l/s Q )d(l/s 6 ) = 0.014 sin 9At(H ohs -H x ). (27) 



O.Ssin^A/i/ias/year, (28) 



which is very similar to its same-shell counterpart (25), except for the sin# modulation. 

Moreover, one has to address the main expected source of noise, to wit the intrinsic peculiar velocities of the 
sources. The variation in angular separation for sources at angular diameter distance Da (measured by the observer) 
and peculiar velocity t> pcc can be estimated as [6] 



Ai7 pe 



500 



km 



' Da 
IGpc 



At \ 
10 years / 



/ias. 



(29) 



This velocity field noise is therefore typically smaller than the experimental uncertainty (especially for large distances) 
and again will be averaged out for many sources. The above relation was further investigated in [50], where it was 
proposed to estimate Da via observations of A t 7 pec due not to voids but by our motion with respect to the CMB. 

Finally, two competing effects induce similar dipolar parallaxes: one is due to our own peculiar velocity and the 
other by a change in aberration of the sky due to the acceleration of the Solar System in the the Milky Way. Both 
effects have nonetheless distinct redshift dependence, which can be used to tell all three apart. Figure 8 depicts the 
three dipolar effects; we will come back to this issue in Section VII B. 



2. Numerical derivation 



As suggestive as the above estimates be, they need confirmation from an exact treatment where the full relativistic 
propagation of light rays is taken into account. We will thus consider in what follows the three LTB void models 
introduced in Sec. II B 1, dubbed Models I, II and cGBH [25]. In all three cases the off-center (physical) distance is 
set to 30 Mpc. Assuming no peculiar velocity of the observer, this distance was shown to be compatible with both 
the CMB dipole [25, 38, 48] and supernovae data [42, 43]. 

To compute numerically the cosmic parallax effect in LTB models, an algorithm was laid down in [6]; it can be 
found in Appendix A. Using this algorithm, we plot in Fig. 7 A t "f for three sources at z = 1, for models I and II as 
well as for the cGBH model and the FRW- like estimate. One can see that the results do not depend sensitively on 
the details of the shell transition and that in both cases the FRW-like estimate gives a reasonable idea of the true 
LTB behavior. 

Fig. 8 illustrates the redshift dependence of the cosmic parallax effect for two sources at the same shell (i.e., same 
redshift) but separated in the sky by 90° (which is the average separation between two sources in an all-sky survey): 
one source is located at £ = —45°, the other at £ = +45°. Also plotted are the two major sources of systematic 
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Figure 7. At 7 for two sources at the same shell, at z = 1, for Model I (full lines), Model II (dashed), the cGBH model (red, 
long-dashed lines) and the FRW-like estimate (dotted) . The horizontal axis is the time-average angular position £ of the source 
in the sky, and can be approximated by the inclination angle 9 to good precision. The plotted lines correspond to a separation 
of 90° in the sky between the sources. The off-center distance is assumed to be 30 Mpc. From Ref. [25]. 
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Figure 8. At7 for two sources at the same shell but separated by 90° as a function of redshift assuming a 30 Mpc off-center 
distance. The dark, brown lines correspond to the cosmic parallax in Models I (full lines) and II (dashed); the red long-dashed 
lines to the cGBH model; the light, blue dotted lines represent 1/40 of the aberration-induced signal (see text), which does 
not depend on redshift; the dark dotted lines stand for the parallax induced by our own peculiar velocity (assumed to be 400 
km/s). Since all effects are dipolar, the curves plotted here are proportional to the amplitude of such dipoles. The actual 
amount of noise depend on the angle between the center of the void and the directions of acceleration and peculiar velocity 
of the measuring instrument. Notice that in Model II the cosmic parallax is zero inside the void, which is expected as H\\ is 
constant inside the void in that model [see Fig. 4 and (23)]. The vanishing cosmic parallax in Model II inside the void is a 
well-understood peculiar feature of that model. From Ref. [25]. 




noise, which will be discussed in Section VII B: our own peculiar velocity and the change in the aberration due to the 
acceleration of the observer. As will be shown, all the effects we are considering are dipolar and the lines in Fig. 8 
are proportional to the amplitudes of such dipoles. Note that both systematics have different z-dependence than the 
cosmic parallax produce in void models, and in principle all three effects can be separated. 



B. Cosmic parallax in homogeneous and anisotropic models 

The Bianchi solutions describing the anisotropic line element were treated as small perturbations to a FRW back- 
ground and its effect on the CMB pattern was studied by [51-56]. These anisotropic models face two main drawbacks: 
in order to fit large scale CMB patterns they sometimes require unrealistic choice of the cosmological parameters; the 
early time inflationary phase isotropises the universe very efficiently, leaving a hardly detectable residual anisotropy. 
However, all these difficulties vanish if the anisotropic expansion is generated only at late time, excited by a domi- 
nating anisotropically stressed dark energy component. In [47] a general treatment of the cosmic parallax in Bianchi 
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Figure 9. Cosmic parallax in Bianchi I models for h x = 0.71, h y = 0.725, h z = 0.72, i.e. E x = -0.012 and E y = 0.009). The 
time interval is At = lOyrs as: (Left) a function of 6 for (f> — A<f> = and A9 = 90°, which corresponds to the (X,Z) plane; and 
(Right) a function of for 6 = 90°, AO = and A0 = 90° which corresponds to the (X,Y) plane. From Ref. [47]. 



I models has been derived, and subsequently connected to simple phenomenological anisotropic dark energy model 
[57]. 

1. Cosmic parallax in Bianchi I models 

The unperturbed metric in Bianchi I models can be written in Cartesian coordinates as: 

ds 2 = -dt 2 + a 2 (t)dx 2 + b 2 (t)dy 2 + c 2 (t)dz 2 , (30) 

where the three expansion rates are defined as Hx = a/a, Hy = b/b and Hz — c/c. Here the derivatives are 
taken with respect to the coordinate time. Bianchi I models exhibit no overall vorticity but shear components 
^xyz — Hx.y.z/H — 1, where H is an effective expansion rate, H = A/A, with A = (abc) 1 ^ 3 . In homogeneous 
and anisotropic models like Bianchi I models we expect a different signal with respect to the dipolar LTB one 
derived in Sec. IIIA1, being less contaminated by systematics (e.g. observer's velocity and acceleration) that could 
mimic the cosmological signal and hence be even more predictive. Let us consider two sources A and B in the 
sky located at physical distance from us observers 0[a.b] — (-^j ^> ^)[a,b] = (-R sin # cos <^>, R sin sin <f>,R cos 9)\a,b] > 
where R = ^JX 2 + Y 2 + Z 2 and (0, 0) are spherical angular coordinates. Their angular separation on the celestial 
sphere reads 

OA • OB = cos 7 = cos #a cos (9b + sin 9 a sin 0g cos A0, (31) 

with A<p = (4>a — 4>b) -If during evolution of the Universe At 7 is different from zero, then a cosmic parallax arises. If 
the expansion is homogeneous but anisotropic, the angular separation between two points, in a At interval, changes 
as: 

— sin7Af7 = sin #,4 cos8B(A t 0B cos A(j> — A 4 6U) + cos 9a sin9B(A t 9A cos A(j> — A t 9s) (32) 
+ sin6U sin6>B sin A(j)(A t (j)B — A t ^). 

Looking at the equation above, some intuitive solutions can be derived. For instance, in the limit of A t (f>A — ^■t't'B — 
4>a = <Pb — the motion is confined on the (X,Z) plane and the cosmic parallax reduces to (Af^ — A t 9s) (see left 
panel in Fig. 9). Similarly, on the (X,Y) plane the signal is (A t 4>A — At4>B) ( see right panel Fig. 9). In figure 9 the 
shear parameters at present are allowed to appreciably deviate from 0. This explains why the cosmic parallax is a 
few orders of magnitude larger than the one in [46] . The main motivation for this will be presented later on. 

In general, the signal will be given by the combination of the anisotropic expansion of the sources and the change 
in curvature induced by the shear on the photon path from the emission to the observer. As for the LTB models, the 
photon trajectory in a Bianchi I Universe will not be radial. However, while in LTB models this effect on the cosmic 
parallax is enhanced by inhomogeneity (although a FRW description of null geodesic has been shown to give a fairly 
good approximation, see Fig 7), in Bianchi I models the geodesic bending for a single source has been shown to amount 
at most to about 7% [47], which allows to adopt the straight geodesies approximation. With this approximation and 
considering the relations (f> = arctan(F/A) and 9 = arccos (Z/^/X 2 + Y 2 + Z 2 ) one can write down the temporal 
evolution equations for the angular coordinates as: 
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Figure 10. Mollweide contour plot for cosmic parallax in Bianchi I models for one source fixed at two different location in 
the sky. Upper panels show the signal for ellipsoidal models (hoz = 0.72 and hox = hoy = 0.71), while in lower panels 
hoy = 0.725. Lighter colours correspond to higher signal and on the horizontal and vertical axes angular coordinates vary in 
the range <j> : [0, 2n] and 8 : [0, 7r], respectively. The time interval is At = lOyrs. From Ref. [47]. 
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(33) 
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where the shear components at present satisfy the transverse condition £ x + S y + T, z = 0. 

Equations (33)-(34) show the quadrupolar behaviour of the cosmic parallax for Bianchi I models in the <j> and 
9 coordinate, respectively. This functional form recovers the one expected for the first non- vanishing multipole 
expansion of the CMB large scale relative temperature anisotropics in Bianchi I model [53]. In general, the pattern 
on the sky of the cosmic parallax signal will be given by plugging the above equations into (32), to wit At 7 = 
At7(#, <fi, A6, A(p, "Eox, Soy, Ho, At), where the only further conjecture is that H does not appreciably vary in At 
[47]. At first order, this seems reasonable for the time intervals under consideration. The signal is independent on the 
redshift, which means for instance that source pairs along the same line of sight undergo the same temporal change 
in their angular separation or that aligned quasars would stay aligned. 

In Fig. 10 a Mollweide projection of the isocontours of the cosmic parallax for peculiar values of the angles 9 and 
<j) with the same choice of parameters of Fig. 9 is shown. In the first case, one of the sources is fixed at the north 
pole and, in the other, one the sources lives in the plane (X,Y). As expected, when the source is at an equatorial 
position the symmetry with respect to the (X,Y) plane is preserved, while when the source is at the north pole a 
symmetry with respect to the (X,Z) plane a cosmic parallax emerges. In a FRW universe the components of the shear 
simultaneously vanish and so does the cosmic parallax. 



2. Cosmic parallax induced by dark energy 

CMB quadrupole has been used to put strong constraints on Bianchi models [53, 55, 56]. In a LCDM Universe, the 
anisotropy parameters scale as the inverse of comoving volume leading so to a natural isotropization of the expansion 
from the recombination up to present with typical limits of the shear parameters of the order ~ 10~ 9 -j- 10 -10 
(resulting in a cosmic parallax signal of order 10~ 4 /ias). However, this constraints must be relaxed in the cases where 
the anisotropic expansion takes place after decoupling, due to, for instance, vector fields representing anisotropic dark 
energy [57]. In [47] the authors analysed the cosmic parallax applied to a specific anisotropic phenomcnological dark 
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energy model in the framework of Bianchi I models [57, 58] (we refer to these papers for details). The anisotropic 
expansion is caused by the anisotropically stressed dark energy fluid whenever its energy density contributes to the 
global energy budget. 

Let us consider a physical model where the Universe expansion is driven by the anisotropically stressed dark energy 
fluid. After recombination, the energy momentum-stress tensor is dominated by the dark matter and dark energy 
components and can be written as: 



T { m )u = diag(-l,w ra ,io m ,w m )p m (35) 
T (dev = diag(-l, w, w + 38, w + 3j)pbe, (36) 

where w m and w are the equation of state parameters of matter and dark energy and the skewness parameters 8 and 
7 can be interpreted as the difference of pressure along the x and y and z axis. Note that the energy-momentum 
tensor (35) is the most general one compatible with the metric (30) [57] . The reason why limits from cosmic parallax 
could be more sensitive with respect to the ones coming from CMB [46] is that the parameters 8 and 7 are allowed 
to grow after the decoupling. Assuming constant anisotropy parameters, an experimental constraint 8 = —0.1 is not 
completely excluded by supernovae data, since it lies on the 2a contours of the 7 — 8 plane, if a prior on w and Q m 
is assumed [57]. More phantom equation of state parameters and/or larger matter densities allow for larger value of 
delta. In addition, and more generally, time dependent 8 and 7 functions, mimicking for example specific minimally 
coupled vector field with double power law potential, can escape these constraints. It can be shown [47] that the 
dynamical solutions for the quantities of interest can be found by expanding around the critical points the generalized 
Friedman equations and the continuity equations for matter and dark energy [57, 58]: 

U' =U(U - 1)[ 7 (3 + R-2S) + 8(3 - 2R + S) + 3(w 
S' =i(9 - R 2 + RS - S 2 ){S[U{8 + 7 + w - w m ) + w, 

R' =i(9 - R 2 + RS - S 2 ){R[U{8 + 7 + w - w m ) + w 

where U = pde/(pde + Pm) and the derivatives are taken with respect to \og(A)/3. In the above equations we have 
introduced: 

R = (a/a-b/b)/H = S x - S F 
S= {a/a-c/c)/H = 2T, x + ^y, 

that naturally define the degree of anisotropy. Since at present the dark energy contribution to the total density is 
about 74% , among the several solutions of the linear system (37) beside the Einstein-de Sitter case 
0), one is interested in the fixed points where the contribution of the dark energy is either dominant: 



- w m )} 

re -l]-6 7 C/} (37) 
n -l]-68U}, 



R. = — -J 5, = — U. = 1, (39) 

d + 7 + u> — 1 0+7 + w — 1 



or in its scaling critical stage: 



3£(<5 + 7 + w) = 3~/(S + 7 + w) = ^ + 7 + ^ U()] 

2(8 2 - 87 + 7 2 ) ' * 2(8 2 - <5 7 + 7 2 ) ' * w 2 - 3(7 - 8) 2 + 2w( 7 + 8) ' 1 ' 

where pde/ Pm = const., i.e., the fractional dark energy contribution to the total energy density is constant. In this 
case, one needs to ensure that the scaling regime in not too much extended in the past, in order to avoid a too long 
accelerated epoch overwhelming the structure formation era [47]. 



3. Forecastings 

Detecting a cosmic parallax for a Bianchi I Universe requires an astrometric instrument with the highest performance 
in detecting quasar positions like the next generation astrometric experiment Gaia (described in Section VII B) . In 
this section we present forecastings using the Fisher matrix formalism and assuming the instrumental specifications of 
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Experiment 


A s 


aacc 


ixt 


(7£ 




Llo detection (oU IVLpcJ 


Gaia 


500,000 


90 /ias 


2.5yrs 


1.5- 10" 3 


190 - 310 Mpc 


< 0.001a 


Gaia+ 


500,000 


50 /xas 


5yrs 


8.3- 10~ 4 


46 - 70 Mpc 


0.03a - 0.24a 


Gaia++ 


1,000,000 


5 fJ,&S 


lOyrs 


6- 10" 5 


1.6 - 2.5 Mpc 


32cr - 49a 



Table III. Specifications adopted for Gaia-like and GaiaH — H experiments, where N s is the total number of sources, a acc is the 
experimental astrometric accuracy, At is the time interval between two measurements, a^ is the la error bar on the present 
shear in the Bianchi I models here considered and the last 2 columns give estimates on detection of different LTB dark energy 
models (we used Models I, II and cGBH, discussed in Section III A 2). The first LTB column gives the la limit on the off-center 
distance X ^ s if one assumes a FRW metric; the second LTB column shows with how many a one would detect a LTB model 
assuming in all cases X b s = 30 Mpc. Model I is the easiest to detect or rule out, Model II the hardest. Note that the error cteo 
does not depend directly on the interval At as it has been factored out in (33)-(34); instead, it depends only on the average 
positional accuracy a acc and the number of sources N s . 



Gaia and 2 other enhanced Gaia-like missions (dubbed Gaia+ and Gaia++) with an average time difference between 
observations At ranging between 2.5 and 10 years [25, 47] and different positional accuracy. Table III lists the details 
of all 3 experiments (we refer to Section VII B for more accuracy about the observational strategy) . 

Since the average time separation between 2 observations is half the mission duration, Gaia (with nominal 5 year 
mission duration) should provide an average At = 2.5 years. Moreover, since accuracy depend strongly on quasar 
magnitude, it is important to estimate how will Gaia's quasar catalogue be distributed in magnitude. Using SDSS 
catalogue as a baseline, a detailed estimate of what the Gaia mission as planned will provide was carried out in Ref. [25], 
and it was found that the average accuracy would be a acc — 90 /ias (see [59] and Section VII B for details) . We also 
include here estimates for two enhanced Gaia-like missions, following Ref. [47]: one with 500,000 quasars uniformly 
distributed on the sphere with a constant average accuracy of a acc = 50 /Ltas and with average time separation between 
2 observations of 5 years (here dubbed Gaia+) and another one, with 1 million quasars, average accuracy a acc = 5 fias 
and average At of 10 years (here dubbed GaiaH — h). Furthermore using average accuracy numbers, it is trivial to 
re-scale the final errors to a different accuracy of any Gaia-like misiion. 

In general Bianchi I models, the cosmic parallax signal depends on four parameters: the average Hubble function 
at present, the time span and the two Hubble normalized anisotropy parameters at present. However, for the allowed 
range of values, contours in the (£ojc,£oy) f rame do not depend on the value of Hq. Moreover, due to the linearity 
of our equations, stretching the time interval between the two measurements or improving the instrumental accuracy 
would result in a trivial scaling on the final constraints presented here. 
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Figure 11. Fisher contours of Cosmic parallax for Gaia+ and GaiaH — H specifications (dashed and solid lines, respectively). The 
double contours identify la and 2cr regions for At = lOyrs. From Ref. [47]. 
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Figure 12. Projected Fisher contours for the skewness dark energy parameters for Gaia+ (upper panels) and GaiaH — h speci- 
fications (lower panels). The double contours identify la and 2a regions for At = lOyrs. The dashed lines represent the case 
of an ellipsoidal universe with w = —1, Uq = 0.74 and 8 = —0.1 (Rq — 0.2) approaching the dark energy dominated critical 
point (where U — 1 and R, ~ 0.3), while the solid lines represent an ellipsoidal universe that has just entered the scaling 
regime, with w = -1, U = 0.74 and 5 = 0.5 (R* ~ -0.5). From Ref. [47]. 



The Fisher matrix is defined as: 

" ^ ol cc <»:,., ' 



where all separations are taken with respect to a reference source and index I runs from 2 up to the number of quasars 
N s to take into account the spherical distances to all other sources. In fact, one should notice that the Gaia accuracy 
positional errors are obtained having already averaged over 2N S coordinates. The Fisher error ellipses are shown in 
Fig 11; the constraints turn out to be of the same order of magnitude of the CMB limits on the shear at decoupling. 
The la errors on T, x and S y turn out to be 8.3 • 10~ 4 and 6 • 10 -5 for Gaia+ and Gaia++, respectively. The 
calculation was performed assuming the null hypothesis (i.e. a Friedmann Robertson Walker isotropic expansion). 
However, due to the linearity of the cosmic parallax equations, the effect of considering alternative hypothesis is 
merely the shift of the ellipses' centers to the values of the new fiducial model. 1 The constraints on the skewness dark 
energy parameters of equations (35)-(36) can be provided by mapping out the Fisher matrix into the new parameter 
space p = (<5,7) via F' — A T FA, where Aij — dJ^oi/dpj [47]. For the scaling solution (40) the error contours are 
shown in the right panels of Fig. 12, for both Gaia+ and Gaia++ configurations. Conversely, if the expansion is 
driven towards a future dark energy dominated solution, equations (39) do not represent the anisotropy parameters 
at present (see [47]). In order to derive a more appropriate functional form for them, the linearized system (37) was 
solved around solution (39) where log(A), fixed to 0, selects the present time. For this second case, results are shown 
in the left panels of Fig. 12. 

The derived constraints are of the order 10 -3 -j- 10~ 4 with a net improvement of about 2 or 3 order of magnitude 
with respect to the limits coming from SNIa data [57]. It is worth noticing that even though the available number 



This is not generally the case with Fisher Matrices, where the area of the ellipses can change appreciably with different fiducial models. 
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of supernovae will drastically increase in the future, it might be hard to improve the constraints at such a high level 
because of the integral dependence of the luminosity distance on the skewness parameters. Therefore the cosmic 
parallax seems to be an ideal candidate for testing the anisotropically stressed dark energy. In this analysis two main 
experimental approximations were assumed. The first one was to consider the experimental covariance matrix as 
diagonal. The off-diagonal correlations is still under investigation by Gaia collaboration and once these are provided 
the Fisher formalism would naturally include them in the final limits. The second approximation was that only 
the statistical errors were considered in our Fisher analysis. We refer to Section. VII B for discussion on possible 
systematic effects such as peculiar velocity of the objects or aberration change induced by our own motion. 

Before concluding we comment briefly on the possibility of testing anisotropy through the accumulated effect on 
distant source (galaxies, quasars, supernovae) distribution. The number density of quasars will change so that the 
number counts should show some level of anisotropy. If sources shifts by O.l^as/year during the dark energy dominated 
regime, then the accumulated shift will be of the order of 1 arcmin in 10 9 years and up to fraction of a degree in the 
time from the beginning of acceleration to now. If the initial distribution is isotropic, this implies that sources in 
one direction will be denser than in a perpendicular direction by roug hly 1/90 « 1CT 2 [47]. This anisotropy might 
be seen as a large-scale feature on the angular correlation function of distant sources, where we expect any intrinsic 
correlation to be negligible. The Poisson noise become negligible for N 3> 10 4 : for instance, a million quasars could be 
sufficient to detect the signal. Although the impact of the selection procedure and galactic extinction is uncertain, this 
crude calculation shows that the real-time effect could be complemented by standard large-scale angular correlation 
methods. 



IV. PECULIAR ACCELERATION 



In Sections II and III we reviewed the time evolution of the cosmological redshift signal (the redshift drift) and of 
the overall angular separations between distant sources (the cosmic parallax), respectively. These are both cosmic 
signatures of a background expansion: the first one is a perfect tool to probe late time acceleration, while the second 
one tracks a possible anisotropic expansion. 

Many spurious effects might spoil these cosmological signals: Earth rotation, relativistic corrections, the acceleration 
of the Sun in the Galaxy and peculiar motion of the source. Of particular interest is the redshift variation induced by 
linear perturbations, which - in analogy with the peculiar velocities arising from the same physical mechanism - we 
dub peculiar acceleration. This signal may contain additional information on the matter clustering amplitude. Also 
worth investigating is the acceleration field in bound structures (galaxies or clusters), which can be used as a probe 
of the local gravitational field and can be a measure of the mass inside structures and help distinguishing between 
competing gravity theories [10, 11]. 

To be fair, the practical observability of these effects is still to be assessed in any detail. The main question is to 
select convenient targets to measure the high resolution redshifts. A possible candidate is the 21cm lines of neutral 
hydrogen clouds: future experiments like SKA plan in fact [60] to measure redshifts with an error of the order of 
10~ 6 . Whether this precision, combined with SKA's high angular resolution, and possibly other atomic or molecular 
lines, is enough to reach the detection threshold remains to be seen. 



A. Peculiar redshift drift in linear approximation 



The redshift is by definition the ratio of the frequency of photons measured at the observer position (O) and the 
emission position (E) . We suppose that the emitter has world line E with unit four velocity Uq and the observer has 
world line O with four velocity u a E . If the light ray is an affinely parameterised null geodesic with tangent vector i a , 
then the expression for the redshift can be written as 



l + z = 



(42) 



where is the photon 4-momentum. In [61] the authors expressed the redshift at first order in perturbations 
assuming a Newtonian gauge g a = — a 2 (l + 2(f>)dri 2 + a 2 (l — 2ip)5 i jdx' l dx : ' and neglecting the effect of gravity waves. 
The resulting expression for the redshift is 



(1 + *) 



a(vo) 



1 + [<(> + eV 



o 
He 



(43) 



22 



where e 1 and v i represent the spatial components of l a and u a , respectively, and prime refers to a derivative with 
respect to the conformal time. Eq. 43 is clearly composed by adding up the cosmological redshift to the gravitational 
redshift expressed in this gauge. The variation of conformal time at E is associated to the variation of conformal time 
at O with a relation that takes into account the difference between cosmic and proper time due to the motion of the 
observer and the emitter, and it reads: Sije — [1 + e • (ve — vo)]<5?7o- The final expression for the redshift drift in this 
coordinate system is the following 

z = z(t] ,z) + C(vo,z,e,x ), (44) 

where (xo,?7o) is the location of the observer, the dot refers to derivative with respect to observer proper time and 
the first order redshift drift in metric perturbations and v is 

C(r] ,z,e,x ) = -</>oHvo,z) + (l + z)[e-v- ^]°. (45) 

Notice that the term e • v is the line of sight peculiar acceleration of the source under consideration. The first term in 
Eq. 45 is clearly generated by the gravitational potential at the observer location, while the second involves respectively 
the variation of the Doppler effect due to the relative motion of the observer and the source = (1 + z)[e ■ v]^) as 
well as the equivalent of the integrated Sachs-Wolfe effect of CMB due to temporal variation of ip between E and O 
(Q = (l + *)Mg). 

On sub-horizon scales tp — <p and V 2 — ^H 2 tt m (z)a 2 5, where Q m (a) is the time dependent matter density 
parameter and the matter density contrast S is proportional to the growth factor D(t). Picking up the growing mode 
D + (t), the time evolution of the gravitational potential reads 

j> = Hcf>[f(t)-l}, (46) 

where f(t) — d log D + /d log a. Using this expression one can easily couple the root mean square fluctuations of <j) to 
00, which in turn are related to the matter density fluctuations crs(z) = J ^J^P$(k, z), where P$(k,z) is the matter 
power spectrum. Assuming that / ~ 1 at the time of the emission and t = is the time of observation, the root mean 
square redshift drift fluctuations caused by the integrated Sachs- Wolfe effect is the following [61] 

< Cj > 1/2 - (1 + z)H [f(0) - 1]^(0). (47) 

In cosmological linear theory it is also possible to express the velocity gradients 9(x,t) — diVi/aH in terms of the 
density fluctuations, such that 9(x,t) = —f(t)S(x). Since the local linear peculiar acceleration is Vi = —Hvi — di<j)/a 
and the Fourier components of the velocity density contrast satisfy k 2 Hvi = —f(t)H 2 aki8k, then the total root mean 
square fluctuations on the redshift drift due to the variation of the Doppler is the combination of two terms, at the 
observer and at the emission location respectively 

< CI > 1/2 (*) = [< do > (*)+ < CIe > W] 1/2 (48) 

< <lo > V2 (*) = (! + «) [ffirno - /(0)] Hl&^) 

< (Ie > 1/2 (*) = [^ ro (z) - f(z)]H(z) 2 as(z), 

with a s (z) = f 0P s (k,z)/k 2 . 

Assuming the explicit function fit) for a simple flat ACDM model and the prescription by Bardeen et al. [62] 
to estimate the matter power spectrum (and as(z) = as(0)), in [61] the authors presented the dependencies of the 
peculiar acceleration contribution to the redshift drift fluctuations compared to differences in cosmological signal 
between a ACDM model and a model with dark energy constant equation of state (Fig. 13). In this particular models 
the dominant component is which rises at a percent level up to z = 4. In contrast, they find that the Sachs- Wolfe 
redshift variance is much smaller, of order < £? > 1 / 2 ~ (1 + z)Hq ■ 10~ 5 (cfr Eq. 47). This means that the ratio of 

and G to the cosmological signal is about 10~ 2 and 10~ 5 , respectively. Obviously, this might not hold for other 
cosmological models. 



B. Peculiar redshift drift in non-linear structures 



At linear level the peculiar acceleration signal acts as a noise over the cosmological signal and one needs a large 
number of sources (and spectral features) to average the latter out. At small distances however the peculiar field 
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Figure 13. Amplitude of the r.m.s. of the systematic errors due to cosmic acceleration effects. The contribution of (i,o 
(dashed line) is subdominant compared to the one of (dotted line). The solid lines represents the difference between a 
standard ACDM model and cosmological models with either w = —0.95 (upper solid line) or w = —0.98 (lower solid line). 
From Ref. [61]. 



dominates over the cosmological one and can be observed directly in objects where the matter density contrast has 
turned non-linear. In two papers [10, 11] it was shown that the peculiar acceleration in nearby clusters and galaxies 
is in fact of the same order of magnitude of the cosmological signal at larger distances and could be measured with 
the same instrumentation. 

Moreover, the peculiar acceleration field can help solving one of the most outstanding issues in astrophysics and 
cosmology, namely the measure of the total mass of clustered structures. The most frequently used methods rely 
on kinematic measurements, where the velocity dispersion of some suitable class of test particles is used to infer the 
virial mass of the object, (e.g. the rotation curves of spiral galaxies and the velocity dispersion of elliptical galaxies). 
This of course requires the assumption of virialization. The same is true for methods based on X-ray emission from 
the hot intra-cluster gas confined by the gravitational potential. On the other hand, as we will show in the next two 
subsections, the measurement of peculiar acceleration acceleration (i.e. the local peculiar redshift drift) is a direct 
probe of the acceleration field and does not assume virialization. 



1. Galaxy clusters 

Let us consider a particle at the spherical-coordinate position (r, a, 7) in a system centered on a cluster. For a 
schematic representation of angle and distances we refer to the left diagram of Fig. 14. The peculiar acceleration of 
a particle reads 

a = V$ = $ r f H r — $ Q a+-<I> 7 7 (49) 

r sin 7 ' r ' 

where the tilded quantities are unit versors and $ is the gravitational potential (Earth's local acceleration is assumed 
to be properly subtracted). The acceleration along the line of sight unit vector s is then 

a s = s • V$. (50) 

For simplicity the potential is assumed to be spherically symmetric, namely $ — $(r). Therefore the line of sight 
acceleration is as follows: 

a s = cos6>'$ r , (51) 

where 9' is as in Fig. 14. In the same figure R c is defined as the cluster distance from observer, r as the particle 
distance from cluster's center and 9 the viewing angle. For small viewing angles 9 (i.e. for r <C R c ) we have that the 
peculiar acceleration along the line of sight is [10] 

a s = cos6»'$ ,. w sin/3$ r | r=fice/cos/ j (52) 
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Figure 14. Definition of geometric quantities used in the text. Left: a schematic diagram of a galaxy cluster geometry. R c is the 
distance of the cluster center from the observer, while r is the distance of a test particle from the cluster center. Right: Kuzmin 
potential in cylindric coordinates for a disc galaxy seen from an observer at distance R g from the galaxy center (GC). The disc 
lies edge-on on the axis R. The curved contour shows an equipotential curve and 3k represents the acceleration experienced 
by a test particle having galactic latitude j3. Note that the acceleration does not point towards the galaxy center. 



and under the asumption of spherical symmetry it follows that 

a 5 « S m/3^W (53) 

where r = R c 0/ cos/3. The assumption of spherical symmetry simplifies the treatment and allows to use the NFW 
profile for the cluster density but it is by no means crucial to the argument and does not imply any assumption on 
the dynamics, i.e. on the degree of virialization. 

Obviously, in nearby clusters of galaxies the peculiar acceleration signal-to-noise should be higher, if we consider the 
cosmological redshift drift as noise. In a cluster the dark matter halo gives the main contribution to the gravitational 
potential . In Ref. [10] a Navarro-Frenk- White matter density profile is assumed 

p( r ) = Wf^rW' ( 54 ) 

where r s = r v /c sets the transition scale from r~ 3 to r , c is a dimensionless parameter called the concentration 
parameter, p cr = 3Hq /8ttG is the critical density at the redshift of the halo, r v is the virial radius inside which the 
mass density equals A c p cr and S c is the characteristic overdensity for the halo given by 

5 C = C ^, (55) 



where 



C = 



log(l + c)-— - . (56) 
1 + cJ 



In addition, A c is the nonlinear density contrast for a virialized object and enters the expression for the mass 
M v = 4/37rr^A c /9 cr . Its value depends on the cosmological model and assuming a ACDM was set to A c = 102 [63]. 
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Then the mass associated with the radius r is 



M (r) = M V C( log (1 + -)- - 

\ V a 1 



(57) 



and consequently 



$(r), 



GM V „/l°g (1 



-C 



1 



(58) 



Considering a time interval At the velocity shift of the particle test due to the peculiar acceleration along the line of 
sight is Av = a s At and turns out to be 



Aw 



GM„ 



-C At sin 



log(l 
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cm n At M v , 
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(59) 



with r = i? c #/ cos/3. For a typical cluster value of C « 1, it turns out therefore that the typical shift for a galaxy 
cluster is of the order of 1 cm/sec , similar to the cosmological value at z s=s 1 [22] . The maximal value depends in 
general on the separation angle 9 between cluster center and galaxy; for 6 — > we obtain 



Av max = 0.44 



cm At M v 



sec lOy 10 14 M, 



lMpc j 2' 



(60) 



If we can observe Av max as a function of 6 for a given cluster then we may measure a combination of the cluster 
parameters M v ,C,r s and compare this with the expected results from numerical simulations. 

The maximum velocity shift of a cluster is plotted in the right-hand panel of Fig. 15 together with the cosmological 
signal for a ACDM with J7a = 0.7. The redshift at which the two signals are similar depends on the mass of the 
cluster and goes from very small (z ~ 0.02) to ps 0.4. In addition, there is a second redshift, around z ~ 2, at 
which the cosmological signal is negligible. However observations at such a distance could not be reliable, not only 
because of the difficulty of precise measurements at this redshift, but also because the zero-crossing of the function 
is cosmology-dependent. For massive and nearby clusters like Coma, the cluster peculiar acceleration is therefore the 
dominating effect. In principle of course the cosmological signal could also be subtracted from the peculiar one after 
averaging over a large number of galaxies. 
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Figure 15. Left: predictions for Av in cm/sec for Coma (R c = lOOMpc, r v = 2.7Mpc , c = 9.4) for four values of 9 up to 
Smax = 0.027 equally spaced, starting from 9 min = 9 max /A: (top to bottom). This value of 9 min corresponds to a radius of 
0.675Mpc. Right: The cosmological velocity shift as a function of redshift for a ACDM model (solid line) and the maximum 
of the velocity shift due to the peculiar acceleration for 9 — 9 max /4: and three different value of the mass of a rich cluster: 
M = 10 14 M Q (short-dashed line), M = 5 • 10 14 M Q (long-dashed line), M = 1O 15 M (long-short-dashed line). From Ref. [10] 



Among the clusters, Coma seems a very suitable candidate: it is the most studied and best known cluster of 
galaxies, it is almost perfectly spherically symmetric, very massive (M v = 1.2 x 10 15 Mq), and close to our Local 
Group (R c = lOOMpc, z = 0.02 ; the cosmological velocity shift is of the order of 10 _1 cm/sec at this redshift [10]). 
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Assuming, as above, a NFW density distribution, the other parameters are: r v — 2.7 Mpc, c = 9.4 and r s = 0.29 Mpc 
[64]. The velocity shift predicted for Coma is shown in Fig. 15. 

In Fig 16 the contour curves corresponding to the same Av max for 9 = 9 max j\ (and different combinations of c 
and the mass, fixing r s to the Coma value) are overplotted on a fitting formula for the concentration parameter as a 
function of the mass derived from N-body simulations [65] . A measurement of Av max at different mass scales would 
provide then a test of the concentration parameter fitting formula [10], i.e. of the relation between C and M v . 

Another observable is the density distribution of galaxies along a line of sight 9 

N(Rc, 9, s) = 2vrp„ ° A . — ^dAv, 61 

v ' 9 \cos/3(Av)J cos 2 (3{Av) dAv y ' 

where p g is the radial distribution of galaxies (possibly derived by the projected number density) . Comparing N with 
the observed numbers one could reconstruct p{r) [10]. 
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Figure 16. Contour plot for the velocity shift of a cluster at 9 max /A as a function of the concentration c and the mass M. 
Higher values correspond to darker regions and the contours levels are (left to right) Aw = (0.1, 0.5, 2, 4, 6, 8, 10, 15)cm/s. The 
short-dashed line is the fitting formula by [65]. From Ref. [10]. 



2. Galaxy 

A similar analysis can be also applied to our own galaxy [10] . Here the effect of the peculiar acceleration is even 
stronger than for distant clusters and, following Ref. [11], we can ask ourselves whether we can use the Galaxy redshift 
drift as a test of gravity. Different gravity theories, in fact, produce a different acceleration field and this could in 
principle be compared to observations. In Ref. [11] the authors modelled the galaxy following two competing gravity 
theories: a baryon disc embedded in a dark matter halo (we call this model disk+CDM halo) and the same disc 
subject to a different Poisson equation beyond a certain scale (the MOND model [66], to be denoted as disk+MOND). 
The question is whether these two completely different descriptions could be distinguished by the velocity shift as an 
observable, given the same rotation curves. Of course the best target to do so is our own Galaxy, where the known 
globular clusters in the stellar halo might be good test particles. The bulge component will be neglected in both 
scenarios, since its spherical symmetry allows one to treat it simply as an additional contribution to the total mass 
for scales outside the bulge. In this section for a schematic representation of angles and distances in the galaxy we 
refer to the left diagram of Fig. 14. The velocity shift is in both scenarios Av — a s ■ At, where a s is the total line of 
sight acceleration caused by disk+CDM halo and disk+MOND, respectively. 

a. Disk+CDM halo If a test particle (e.g. stars, gas) orbits on the disc, its peculiar acceleration is affected only 
by the mass embedded within its distance from the centre, due to the symmetry of the distribution, and it is the same 
as if this mass was totally concentrated in the centre. 
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The disc component can be modeled as the so-called Kuzmin disc, namely a disc with superficial density 

w \ hM 
S(r) 



2it{R 2 + h 2 ) 3 / 2 ' 



(62) 



where M is the total disc mass, R and z are cylindrical coordinates and h is the scale length of the disc. The 
two-parameter Newtonian gravitational potential outside the disc is 
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[ff2 + Q z \+ h )2] 1/2 < 



(63) 



and the equipotential surfaces are concentric spheres centered at ±/i (see left diagram of Fig. 14). The test particle 
acceleration then reads 
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(64) 



with \z\ accounting for the acceleration field above and below the disc. The force field of a Kuzmin disc no longer 
converges towards the origin of axis. Consequently, the projection angle 9' must be corrected by 7, namely the angle 
between the acceleration and the radial direction r. On the other hand the peculiar acceleration of a test particle 
outside the disc, e.g. a globular cluster, is not only affected by the gravitational potential generated by the disc, 
but also by the possible presence of a dark matter halo needed to explain the observed rotation curves. In standard 
Newtonian mechanics the halo potential is modelled as a logarithmic function 
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(65) 



where R c is the scalelength, q is the halo flattening (q — 1 recovers spherical symmetry) and vq is the asymptotic 
value of the velocity at large radii. 

The Kuzmin and the halo accelerations must be projected along the line of sight and then added together. While 
for spherical symmetric logarithmic potential the acceleration is radial and the angle between the line of sight and 
r is simply 9' (as in Fig. 14), as mentioned before, the projection angle for the Kuzmin acceleration does not point 
towards the origin. One clearly has cos 9' sin (3. Using r cos (3 ~ R g 9, z ~ ±V~r 
sight accelerations read [10] 
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b. Disk+MOND The Modified Newtonian Dynamics (MOND) paradigm was first proposed by Milgrom to rec- 
oncile discrepancies between general relativity and galaxy scale dynamics [66]. As presented in its first version, 
MOND by itself violates conservation of momentum and energy. The subsequent Bekenstein-Milgrom formulation of 
MOND [67] leaves the Newtonian law of motion intact and modifies the standard Poisson equation for the Newtonian 
gravitational potential as follows: 
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AnGp, and ao is a scale that was estimated by Milgrom to 



where ip is the MOND gravitational potential, V 2 
be a = 1.2 • 10 
fj,(\a\/a ) = 1 for 

Eq. 68 is a non-linear equation which might be hard to solve analytically, except for a class of symmetric configura- 
tions. In [68], a class of disc-galaxy models has been introduced for which exact solutions of the MOND field equation 
exist. Subtracting the usual Poisson equation from the MOND equation (68) one has 



10 m/s 2 . The original shape of /1 that helps rendering the right profile of rotational velocities is 
o| » |oo| and fj,(\a\/a ) = \a\/a for \a\ <C |a |. 
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For configurations with spherical, cylindrical or plane symmetry the relation between the MOND field and the New- 
tonian field becomes 



V a n / 



JV- 



(70) 



This equation permits a straightforward relation between the two potentials, such that by assuming a matter density 
distribution, one can solve the Poisson equation for the Newtonian potential and then invert it to get the MOND one. 

The function n defined in MOND has been reformulated during the years, moving from a step function as the 
aforementioned type to, e.g., fi(x) = x/\l + x 2 . As pointed out in [68], the function I(x) — Xji{x) is therefore 
monotonic and invertible, and the inverse is related to the function v(y) = I~ 1 (y)/y- Eq. 70 now reads as 
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The exact MOND solution is then given by Eq. (71), with 

a-M = - yip = a a I — — 

where clm and ajy are the MOND and Newtonian acceleration, respectively [11]. 
For fx(x) = x/y/l + x 2 , then v{y) = [1/2 + y/y~ 2 + 1/4] and consequently 
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The expression for the peculiar acceleration in MOND, outside the disc, is as follows: 

a s ,M = a M - sin (/9=F7), 

with 
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The different shape of the peculiar acceleration between Newtonian and MOND configurations are shown in Fig. 17. 
The predicted velocity shift reaches a maximum value along the line of sight, Av max , that depends on the gravitational 
configuration. The maximum also depends on the distance from the galactic centre, increasing with decreasing 
distance. By estimating Av ma x from a sample of test particles (e.g. globular clusters) outside the disc of the galaxy 
along each line of sight, one might therefore hope to distinguish the Newtonian and the MOND gravitational potentials. 
In our Galaxy, knowing precisely the coordinates of test particles one might trace out the distribution of the signal. 

In this case we clearly are internal observers. The projected acceleration of the Sun on the disc, which has its 
maximun value on the order of a few cm/s, is assumed to be substracted, and the cluster-centred acceleration of 
single stars in globular clusters is assumed to be averaged out. The time scales are short enough that the motion of 
the Milky way within the Local Group and the Local Group relative to Virgo, the Great Attractor and other distant 
cosmic structures can be considered approximatively constant. In [11] the velocity shift signal has been forecasted for 
150 globular clusters catalogue in our Milky way. 

The effect of being off-centred observers combined with the non-radial direction of the acceleration in Kuzmin disc 
makes the pattern of the effect across the Galaxy non trivial. However, while again the signal reaches the maximum 
value close to the Galactic centre for both scenarios, the difference signal between the two Galaxy models is stronger at 
high Galactic longitude (close to it) where the spherical CDM halo is more important. This can be better seen by eye 
in Fig. 18, where the distribution of the peculiar acceleration in Sun-centred cartesian coordinates has been mapped 
out . Unfortunately, the maximum difference is only of order 1 cm/s, making it quite problematic to observe. However, 
supposing that such an accuracy will be achieved in the future, it is interesting to note that the two different scenarios 
result in a distinct morphology of the signal distribution, that could be used as a signature to identify the actual 
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Figure 17. Left: velocity shift curves caused by peculiar accelerations as a function of ft for Newtonian configuration (thin 
lines, CDM halo+Kuzmin disc) and MOND (thick lines). For both configurations curves are drawn for three different values 
of 9, from top to bottom: {0.05,0.06,0.07}, corresponding to distances from the centre {40,48,56} kpc. -.Right: same as the 
left panel, but for three different values of 6, from top to bottom: {0.01, 0.03, 0.05}, corresponding to distances from the centre 
{8,24,40} kpc. From Ref. [11]. 



potential (compare Ref. [11] for a Mollweide projection of the signals in the sky). Here the authors also provided 
a top-ten list of the globular clusters exhibiting the highest velocity shift signal and the first ten with the highest 
difference signal between the CDM halo configuration and MOND. The highest signal among the globular clusters 
predicted for MOND configuration, CDM halo and their difference is 21.5cm/s, 20.8cm/s and 1.18cm/s, respectively 
(in 15 yrs). 

The advantage of adopting globular cluster as test particles is related to the fact that they not only are bright, but 
more importantly composed of thousands of stars and their position is fairly well known; nevertheless, in principle one 
might consider using other probes, like high velocity clouds that consist of neutral hydrogen high velocity interstellar 
regions or HI regions. Moreover, one could also use the acceleration in the pulsar timing as another probe of the 
acceleration field in the Galaxy (e.g. by measuring the Doppler shift in the pulsar timing). The practical feasibility 
of such probes is however still to be demonstrated. 



V. PROPER ACCELERATION 

As we mentioned in the Introduction, the fourth application of real-time cosmology is to the velocity variations 
in the direction perpendicular to the line of sight (i.e. transverse). When this effect is observed in test particles in 
bound systems, for example within our own Galaxy, we call it proper acceleration (since it represents a variation in 
the proper motion of the source) to distinguish it from the cosmic parallax. Let us stress that the material presented 
in this section is set forth here for the first time and it is somewhat preliminary. 

Let us assume that a body at distance r s from the Sun has transverse velocity Vtx- Then during the time AtM, in 
which an astrometric mission like GAIA is flying, its angular position will change by s\ — VxiAtM /r s - Then another 
GAIA-like mission At$ years later will see the change s 2 = v^At^j l r s (assuming the same duration for simplicity). 
If vt2 7^ vtx then we will see a proper acceleration. If r s changes also during the interval Ats, we see also an apparent 
proper acceleration (see below). 

If the body has a (slowly varying) transverse acceleration we have vti ~ v ti + cixAtg and therefore the variation 
in proper motion will be 

A t s = s 2 - si = a T At M At s /r s . (77) 

We can produce a quick estimate of the effect assuming radial symmetry around the galactic center. If the body 
lies in the direction r s , is at distance r from the galactic center, and if the galactic center is in direction r c , then 
r 2 = r\ + r 2 s — 2r c r s cos a, a being the angle between r c and r s (see Fig. 19). Then we may put 

y2 ■ a' v2rc ■ t>7Q\ 
fly = — sin ft = — -sma (78) 
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Figure 18. The velocity shift signal (for At = 15 years) expected in the MOND configuration (top left), in the newtonian 
configuration (top right) and for their difference (bottom) for the Milky Way globular clusters taken from Harris (1996) [69], 
plotted in cartesian coordinates (kpc units). From Ref. [11]. 



9' being the angle between the line of sight and f . Then we have 

v 2 s\naAt M At s r c 

A * s = 2 i 2 o ■ ( 79 ) 

r c + r s ~ iTc^s cos a r s 

For a star at r s = lOOpc, with v = 300Km/sec and a — ir/2, and adopting Atj^ = 5yrs and At$ — lOyrs, we 
obtain A t s ~ 1/xas. Now suppose GAIA observes N objects with the same acceleration (e.g. within a cluster or in a 
small region of space) with precision P/\fN{ias each. Therefore it will be able to observe the proper acceleration if 
A t s > P^as/%/]V or equivalently 

r s \ /300km/sec\ 2 / 10yr\ / 5yr 
N = QMP< w 1 ' ' J u * 



/ rs\ /dUUkm/secy /lUyrx / 5yr \ 

lioo^Jl v J Va^JIa^J (80) 

x 1 + (r s /r c ) 2 - 2(r s /r e )cosa 
sin a 

where we assumed r c — 8.33kpc. It appears therefore that for bodies in the disk near the Sun (e.g. an open cluster 
at lOOpc) and a sa 7r/2, the effect could be visible if P ~ 10/xas and one has ~ 10 2 stars to average, which seems 
not impossible. To reach a distance of cikpc with astrometry precision Pfias one would need to average over roughly 
lOOcPP 2 stars. In order to realize the measurement with a single GAIA mission one should put Aij\/ 2yr and 
Ats ~ 5y r - Then one gets N ~ 2500d 2 P 2 which appears quite feasible. 

The measurement will be subject to several source of errors. We consider here in particular the apparent acceleration 
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Figure 19. Kuzmin disc potential in cartesian coordinates in the Milky Way. The disc lies in the plane (x,y). The star 
(S) represents the Sun and angles are labelled with a hat. Curved contours show two equipotential levels; 3k represent the 
acceleration of a test particle at galactic latitude 0. 



A t s a induced by a change in r s . We have 

a - a , / I 1 s A , Ar s At M At s 

A t s a = v T iAt M { ) ~ v T iAt M ^r = vtivr. 7. (82) 

r S 2 T s i rj. rj. 

which is r /r s (vRVTi/v\ + uf-x) ti mes bigger than the true acceleration. If vr ~ vti and the object is w lOOpc from 
the Sun then A t s a / A t s ~ (r c /r s ) 100). A signal-to-noise of order unity is reached for r s r Cl but then we would 
need roughly 10 stars to average in order to detect the effect. However the apparent acceleration can be measured 
with good precision by estimating vr with Doppler measurements and r s via parallax distance. The relative error on 
the estimate of A t s a is given by summing the variance on the redshift and the variance on the distance: 



er„ — u„ + ov 



(83) 



GAIA will measure radial velocities with an error of 15 km/sec, which amounts to a z « 10 _1 , and parallax distance 
with an error smaller than 10%, so again a r ^ ps 10 _1 . Therefore we can estimate a a w 0.1. Then we need to have 
a ratio A t s a /A t s < 10 to ensure detection of the true proper acceleration. This gives an lower limit of r s > lkpc 
roughly. Therefore we need roughly at least 10 4 stars to average over if P = 10/J,as. 

As in Section IV we can now try to see whether proper acceleration can help distinguishing between the acceleration 
field induced by Newtonian and MOND gravity. For a general gravitational potential the signal is 

. _ $ r sin a At M At s 

tS ~ r s [l + (r s /r c )2_ 2 (r s /r c )cos«] 1 /2- ^ > 

Adopting the configuration introduced in Sec. IV B 2 for a disk plus a CDM halo depicted in Fig. 19 the transverse 
signal we are after is then: 

A t s = (a ttK + a t . L )At M At s /r s . (85) 

The Kuzmin and the halo accelerations must be projected on the transverse plane and then added together. While 
for spherical symmetric logarithmic potential the acceleration is radial and the angle between the line of sight and 
r is simply 9' ' , the projection angle for the Kuzmin acceleration, again, does not point towards the origin. The two 
transverse accelerations read 

MG 

at ' K= [R* + (\z\ + h)>] Sln(g,±7) (86) 

I 4. 5i 

a t l = v 2 a — i ^-5- sin 9'. (87) 
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All the parameters of the gravitational potential will be fixed to the values of [70]: M = 1.2 • 10 11 M Q , h = 4.5 kpc, 
vq = 175 km/s, R c = 13 kpc and q = 1. 

Let us then write the contribution from proper acceleration to velocity shift in terms of galactic coordinates, in 
particular r s , the galactic latitude b and longitude I. First, we have to move from cylindrical coordinates (z,R,ip) to 
galactic coordinates (r s , b, I), and neglecting ip due to the axisymmetry of the potential, the transformation reads 

z = r s sin b; R = \J r| + r 2 cos 2 b — 2r c r s cos b cos I , (88) 

The line of sight corresponds to the direction of the distance from the object to us, i.e. r s . The projection angle 9' 
is related to the new coordinate system as follows: 

cos9' = ^ c -, (89) 

2r s r 

where r 2 = z 2 + R 2 . In addition, as already mentioned, for the Kuzmin potential the correction to the non-radial 
direction of the acceleration must again be taken into account, so that the projection angle changes to cos (9' ±7), 
and 

t 2 ~\~ x 2 — h 2 

7 = arccos ; x 2 = r 2 + h 2 - 2rhsin(3. (90) 

2rx 

with 

/3 = arccos c 1. (91) 

2rr c 

Using these relations we can find the locations where the signal is strongest. We find that the maximal newtonian signal 
is 0.003[A£jv/Ais/(5yrs 2 )]/ias at r s — 3kpc, and decreases to 0.0002[A£a/ Ais/(5yrs 2 )]pas at r s = 8kpc. For the differ- 
ence |MOND-newton| at the same distances we find 0.0018[Ai A /A<s/(5yrs 2 )]^as and 0.00006[Atj\/Ai,g/(5yrs 2 )]/xas, 
respectively. 



VI. COSMIC MICROWAVE BACKGROUND RADIATION 



Recently, several authors also entertained the possibility of detecting real-time variations of the CMB angular power 
spectrum [13, 14, 71]. The future evolution of the temperature anisotropy power spectrum should be affected by three 
effects: i) the amplitude should reduce due to the \/a scaling of the mean temperature; ii) the peaks and the minimum 
should shift to smaller scales, because the radius of the last scattering surface will grow and the same comoving size 
will be seen under a smaller viewing angle; iii) in presence of either curvature or dark energy the late integrated 
Sachs- Wolfe (ISW) signal would increase with respect of the first peak. In addition our own velocity with respect to 
the CMB frame evolves, hence generating a temporal change in the dipole. 

As it will be clear later on, the prospect for the detection of this signal are at the moment rather weak. 



A. Monopole 

The first mechanism responsible for the temporal change of the CMB anisotropics is the evolution in time of 
the mean temperature T. Assuming that at present time the CMB radiation energy density is with a fairly good 
approximation decoupled from the rest of the Universe, its conservation equation holds, i.e. p 7 = —AHp~ r Due to 
the Stefan-Boltzmann law the energy flux density is directly proportional to the fourth power of the black body's 
thermodynamic temperature, p 7 oc T 4 , which leads to the evolution equation for the mean temperature: 

t = —HT. (92) 

In [14] it has been shown how, assuming Hq = 73km s _1 Mpc -1 , this temperature will drop by 1/iK in 5000 yr. 

Indeed T will become during t he expansion harder and harder to detect, but there is also a model dependent 
theoretical limit Tu m = \J A/12n 2 below which the CMB will become dominated by the thermal noise of the de Sitter 
background. Considering the scaling law of the radiation energy density p 1 oc (a/ao)~ 4 it turns out that we must 
wait until a/ao ~ 10 30 before this fundamental limit has been reached; in a ACDM model this corresponds roughly 
to t = lTyr [14]. This monopole effect is clearly out of reach. 
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B. Dipole 



The observed dipole in the CMB anisotropy power spectrum is believed to be caused by the Doppler effect arising 
from our peculiar velocity, v with respect to a frame of reference external to the Milky Way for example, the 
International Celestial Reference Frame. As we have mentioned in the previous sections, due to the local acceleration, 
this peculiar velocity, and hence the dipole, is expected to evolve with time. 

Today we can detect the modulation of the Earth's motion around the Sun; in the future, with increasing satellite 
sensitivity, we may be able to observe the Sun's motion around the Galaxy. In [14] the authors assumed that the 
motion of the Sun around the Milky Way is simply a tangential speed of 220 kms -1 at a distance of 8.5 kpc. Using 
the current observed value of v to infer the velocity of the Galactic centre with respect to the CMB rest frame, the 
time dependent Sun-CMB velocity vector then reads 
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(93) 



where the Galactic orbital period is T — 2.35 x 10 8 yr. 

The most straightforward method to work out whether a temporal change in the dipole is detectable or not, is to 
compute a sky map of the dipole at two epochs. Clearly then, if the temperature variance of the difference map is 
greater than the experimental noise , a detection is likely to be expected. The variance of the difference map, denoted 
in [14] by Cs, is Cs = [(8v x ) 2 + (5v y ) 2 + (5v z ) 2 ]/ (Air), where Sv is the difference of the Sun-CMB dipole vector between 
the two observations. Using Eq. (93), and converting to fractional temperature variations, the signal variance of the 
changing dipole is then 
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(94) 



In Section VII C whether or not this might be achieved by future observations will be discussed in more details. 

It should stressed out that the above discussion relies on the stated usual assumption that the dipole is mostly (or 
completely) due to a peculiar velocity. This, however, is yet to be proven observationally and there still remains a 
possibility that the dipole is due (in large or small part) to more exotic origins (see [72, 73] for a brief review of the 
possibilities). A cross-check that the dipole is really due to a peculiar velocity can nevertheless be made by looking 
for correlations between ai m [defined below in (96)] and a^ m (see [72-75]). 



C. Primary anisotropics 

Besides the time evolution of the dipole, the core of the information about CMB is expanded along all the anisotropy 
scales through the whole power spectrum. Therefore in principle it is much less trivial to derive the temporal progress 
of the other multiples and their correlations. However the authors in [14] have shown how under simple assumptions 
it is possible to determine a scaling relation between power spectra at different times. The main approximation is 
that all of the CMB radiation was emitted from the last scattering surface (LSS) at some instant when electrons and 
photons decoupled, and then propagated freely, supposing that LSS has vanishing thickness. 

The temperature fluctuations observed in a direction n of the sky can be expanded in spherical harmonics as 

6J ^lll = J2 a Ur)YUn), (95) 

where r is the conformal time. All the details about the map of the CMB at a time r are encoded in the expansion 
coefficients ae m , that are related to the primordial perturbations 1Z and the linear transfer function T as follows: 

a em (r) = / kdkK em (k)T(k,e,T). (96) 
The statistical properties of the primary anisotropy ae m s are then directly connected to 

(K hn (k)K}, m , (k')) = 2ir 2 5(k - k')^-S w S mm ' , (97) 

where Pu{k) is the Gaussian curvature primordial power spectrum. Their variance Ci(t) as a function of time, defined 
such that < ai m {T)a* im {T) >— Cf,{r)8u' 8mm' j results to be 

C £ (r) ee 4tt / ^-P n (k)T 2 (k,e,r). (98) 
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From the previous equation it is evident that all the temporal information of the CMB primary anisotropy power 
spectrum is embedded in the linear transfer function; as the comoving radius of the LSS increases with time, every 
wavelength will span a smaller angle. 

Since for large £ one has T 2 {k, I', r') ~ (r 2 /r' 2 )T 2 (fc, £, r) ([14] ), where the scale £ has shifted in time to I 1 = £t'/t, 
the final geometrical scaling relation for the power spectrum is 

e ,2 c e ,(T')~e 2 c e (T). (99) 

In particular, the Sachs- Wolfe plateau for a scale invariant power spectrum remains constant in time, as £(£ + l)Cg is 
independent of £ for a Harrison-Zel'dovich initial spectrum. The polarization spectra will scale as well as Eq. (99) , 
since the main contribution is sourced at LSS. However, minor contributions come from reionization, whose comoving 
radius will then affect the polarization scaling expression. In addition, in a ACDM universe observers endure a future 
event horizon, that is a finit Tf and consequently a maximum £f. 

The measurable real time quantity is the power spectrum difference A T Cg = CV(r') — Ci(t) at two different times. 
For small time spans At it reads [14] 

A T C e c ^-C e (r)Ar = - — (l^l + 2Ce (r)) (100) 

and as all the real time observables it is proportional to the time interval. This temporal power spectrum variation 
can have either negative or positive sign and vanishes for the scale-invariant Sachs- Wolf plateau and at an acoustic 
peak. This happens because assuming statistical homogeneity of space the quantity £(£ + l)Ci does not depend on £ 
in the Sachs- Wolfe plateau for a scale-invariant power spectrum, which is valid for the entire acoustic peak structure. 
The fluctuations of the same physical size in the LSS visible at a time r (corresponding to a comoving distance r) on 
angular scale 9 will be visible at a time r' (corresponding to r') at a smaller angular scale 6' = 9r/r' [14]. 

D. ISW 

The integrated Sachs- Wolfe (ISW) contributions to total anisotropy power spectrum comprise an early-ISW gener- 
ated by the exit from the radiation era and a late-ISW caused by the onset of the dark energy component domination 
at late times. The former takes place during a time that is relatively close to the LSS, hence superposing coherently 
to the primary anisotropies; therefore it is supposed to scale as Eq. (99). The large scale temperature fluctuations 
arising from the late-ISW are due to the small redshift varying gravitational potential, hence it can be described by 
an extra transfer function T(k,£,r). The ISW power spectrum can then be written as [14] 

C [™ ir )=4* J f P.( fc )T 7 W( fe ,^) = g^^dr'(^) 2 (r-r'), (101) 

where g(r) is the perturbation growth function. The time derivative of the previous expression gives the temporal 
variation of the late-ISW effect, to be added to Eq. (100): 

so that the total CMB shift is the sum of Eq. (100) and (102). Left panel of Fig. 20 shows the power spectrum at 
different future times [13]. Note the progressive enhancement of the ISW with respect to the acoustic peaks. Of 
course this is completely model dependent, in that in a Universe with no cosmological constant we do not expect the 
large scale plateau to change. On the right panel a numerical derivation of the CMB temporal shift calculated with 
the Boltzmann code CAMB is compared to the analytical expression, which matches the former reasonably well [14]. 



E. CMB maps 

Knowing the Ci(t) it is possible to draw a series of maps at different times; assuming a gaussian random field 
all the information is encoded in the Ci as they are the variance of the ai m s. In [13] the authors have generated 
random complex values for the multipole coefficients with variance Cg. Since the different maps are supposed to be 
correlated in times, the ai m s need to be derived according to the covariance matrix C/ =< ai m (Ti)a^ m (Tj) > and 
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Figure 20. Left the temperature angular power spectrum of the CMB at several representative time steps into the future taken 
from [13]. The scale factors and conformal distances go from (1, 0) to (602, 1-10 1) counting from present time (confront [13]). 
Right: Absolute value of the difference in the CMB power spectra 1(1 + l)\5Ci / (2ir) between a b s = 1 and a' obs = 1 + 8a. The 
solid curves are from numerical integration in [14], and from top to bottom denote 5a = 0.01, 0.001, and 10 -4 . The dashed 
curve was calculated using the analytical expression for 5a = 10 -4 . 
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Figure 21. CMB maps (units on the bar are uK) for a ACDM model at 4 different future times; from top left clockwise: present, 
1 ■ 10 9 yr, 5 ■ 10 9 yr, 40 ■ 10 9 yr in the future. From Ref. [13]. 



can be expressed as ai m = M^^/D^x, where x is a vector of complex random deviates and Dg is the diagonalized 
covariance matrix, such that Cg = M^D^M|. In Fig. 21 four CMB maps at four future time slices are shown [13] for 
a WMAP best fit cosmological parameter set. The evolution first affects the small angular scales, shifting them to 
smaller sizes due to the expansion of the CMB photosphere, then, for larger time span, the larger scales through an 
evident increase of the ISW. 

For small time intervals the maps seems to remain highly correlated, while as the time span increases the correlation 
vanishes first on smaller scales and then also on larger scales. According to [13] for roughly 10 10 yr the future sky 
decorrelates with present time sky. A detailed analysis of these time correlation can also be found in [14]. In 
Section VII C the detectability of this real time observable will be examined. 
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VII. OBSERVATIONS 



In this Section we review some of the experimental requirements needed to perform the previously described real- 
time observations. Actually, this issue has barely started to be properly investigated and is far from being settled. 
All the knowledge to date clearly arises from numerical simulations and forecasts. 

In particular, we would like to point out that in general one might be misled by the meaning of the time span 
At adopted by all the works on real-time cosmology. For example, the time interval that enters Eqs. (5), (22), (59) 
and (100) is the time elapsing between two independent measurements. This should not be confused with the effective 
time needed to perform the observations. In the next subsection the question will be addressed in more details. 



A. The redshift drift and ultra-stable spectrographs 



Recently, two high-precision spectrographs were proposed which could in principle be used for measuring the redshift 
drift: the Cosmic Dynamics Experiment (CODEX) [9, 19, 76] at the European Extremely Large Telescope (E- 
ELT) [77] and the Echelle Spectrograph for PREcision Super Stable Observations (ESPRESSO) [9, 76, 78] at the 
Very Large Telescope array (VLT). Although proposed later, ESPRESSO would serve as a prototype implementation 
on the technology behind CODEX as part of its feasibility studies and could be operational several years before 
CODEX [9, 76]. 

The possibility of detecting the redshift drift with CODEX was analyzed in a number of papers [9, 20, 22, 76]. The 
achievable accuracy on a& v by the CODEX experiment was estimated (through Monte Carlo simulations) [19] to be 

«. - 1-35 (MV Y*9»V Vli^Vws, (103, 



^2370/ V 30 
with 

q = -1.7 for z < 4 , q = -0.9 for z > 4 , (104) 

where S/N is the signal-to-noise ratio per pixel, Nqso is the total number of quasar spectra observed and zqso their 
redshift. Note also that the error pre-factor 1.35 corresponds to using all available absorption lines, including metal 
lines; using only Lyman-a lines enlarges this pre-factor to 2 [9] . 

Achieving a high S /N level (we refer to Appendix B for a more precise formulation of the S /N level) requires long 
integration times i; nt . In [22], an estimate was performed assuming a CODEX-like experiment, coupled to a 42 m 
telescope with approximately 20% total efficiency, could measure 40 spectra with a S/N of 3000 in about 15 years 
(assuming a 20% use of the telescope, and a 90% of actually usable data). In fact, again according to [19], a CODEX- 
like experiment, coupled to a 60 m telescope with approximately 20% total efficiency, would give a cumulative S/N 
of 12000 for a single QSO, requiring roughly 125 hours of observation to get a S/N of 3000 on that spectrum. Then, 
starting with 10 hours of observation per night, and taking into account a 20% use of the telescope, and a 90% of 
actually usable data, one finds that 40 spectra can be measured with that S/N in roughly 7.6 years (this time would 
actually increase to about 15 years if the telescope aperture is 42 m instead of 60 m) . This integration time is not 
negligible with respect to the interval At b s over which we expect to measure the redshift drift, as assumed in (103). 
Furthermore, it has been claimed in [9] that in principle it would be preferable to spread the observations more 
evenly over At t, s , inserting an unused temporal window between the two measurements; in addition the same authors 
conclude that the best strategy to minimize the errors would be to concentrate as much as possible the telescope time 
in both the beginning and ending of At (, s . Either way, the error estimate (103) is changed somewhat, but never by 
more than a factor 2. However, estimating such a correction depends on the details of the observational strategy. 

In Ref. [25] a compromise strategy was assumed: a three-period observation, each of At }, s /3 duration, and with 
observations contained in the first and third periods. Doing so means that the effective At b s for the redshift drift 
is 2At (, s /3. It is important to note that a larger observational time- frame allows not only for a larger redshift drift 
(which is linear in time) but also for smaller error bars, as more photons are collected and, therefore, a higher S/N 

3/2 

can be achieved. In other words, the "effective signal" increases with At obs if one assumes a proportional telescope 
time is maintained. 



B. Cosmic Parallax and high accuracy astrometric missions 



A realistic possibility of observing the cosmic parallax is offered by the forthcoming astrometric Gaia mission. Gaia 
will produce in five years a full-sky map of roughly 500, 000 quasars with positional error p between 10 to 200 fias (for 
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quasars with magnitude V — 15 to 20). To compare observations to Gaia we need to evaluate the average A t j with 
At = 5 years and N sources. The final Gaia error p is obtained by fitting 2N independent coordinates from iV 2 /2 
angular separation measures; the average positional error on the entire sky will scale therefore as (27V) _1//2 . Over one 
hemisphere we can therefore estimate that the error scales as p/y/N. Since the average angular separation of random 
points on a sphere is tt/2, the average of A t j can be estimated simply as A t "f(9 = 7r/2). Numerically it was found 
A(7(7r/2) = lOs/xas, with little dependence on Az [6]. Therefore Gaia can see the parallax \ipj\J~N < lOs/ias. For 
s = 4.4 • 10 -3 (i.e. the current CMB limit) and p = 30 /Ltas we need N > 450, 000 sources: this shows that Gaia 
can constrain the cosmic anisotropy to CMB levels. An enhanced Gaia mission with At — 10 years (or two missions 
5 years apart), p = 10 /xas and N = 10 6 would give s < 5 ■ 10~ 4 , i.e. r < 1 Mpc if we assume the sources are at 2 Gpc. 

Two local effects induce spurious parallaxes: one (of the order of 0.1/ias year" 1 ) is induced by our own peculiar 
velocity and the other (of the order of 4-5/iasyear _1 [79] by a changing aberration. The latter is the secular aberration 
drift, that is an alteration of the velocity of distant objects (like quasars) driven by the acceleration of the Solar 
system barycenter in the Milky Way. Both produce a dipolar signal, just like a LTB: however, the peculiar velocity 
parallax decreases monotonically with the angular diameter distance, while the aberration is independent of distance 
and is directed toward the Galactic center [79] . Recently, authors in [80] analysed geodetic and astrometric VLBI data 
from a sample of quasars of 1979-2010 and estimated the acceleration vector, together with a non significant global 
rotation, to have an amplitude of 5.8±1.4 /xas/yr and an almost constant dipolar function of redshift. In contrast, the 
LTB signal has a characteristic non-trivial dependence on redshift: for the models investigated here it is vanishingly 
small inside the void, large near the edge, decreasing at large distances. It is therefore possible in principle to subtract 
the cosmic signal from the local one, for instance estimating the local effects from sources inside the void, including 
Milky Way stars. A detailed calculation needs a careful simulation of experimental settings (including possibly effects 
like source photocenter jitter and relativistic light deflection by solar system bodies) which is outside the scope of this 
review. Moreover, more general anisotropic models will not produce a simple dipole. One example is offered by the 
cosmic parallax induced in Bianchi I models [47], as analysed in Section IIIB. 

C. Cosmic Microwave Background accuracy 

The issue of realistically measuring the temporal evolution of the CMB has been addressed in details in two 
papers [13, 71]. The comparison between the accuracy reached by an experiment and the variation of the CMB can 
involve either the power spectrum or the map. First of all it is important to notice that since we are looking at the 
evolution of a single sky realisation the intrinsic uncertainty of the signal due to the cosmic variance should not be 
taken into account. All that matters is the instrumental accuracy and the possible foreground contamination. In 
addition a portion of information is embedded in the time-time correlation that render the simple power spectrum 
less statistically significant than usual. One can then define the difference between a£ m s at two different times Dg =< 
\a,i m {Ti) — ai m {Tj)\ 2 >— Ci{Ti) + Cg(Tj) — 2C\ J . In [71] the authors simulated a full sky dataset with no noise at present 
time. The predicted evolution of the dipole due to galactic motion is (ai-i, aio, a>n) — 740.5 • 10~ 6 (cos z,0, sins) + 
(—1556.5, 919.3, —87.8) • 10 -6 , with x = 2irt/T, T being the galactic rotation period. They then use this prediction to 

build a second noiseless map at some time in the future. The estimator is x 2 = J2i=r [Ti/T(t) —Ti/To] 2 / (2Cn)-, where 
Cjsi =< (ST{ lmse /T) 2 >= a 2 ix is the noise uniform variance. In the noise dominated regime this statistic reduces to 
the signal-to-noise ratio related to the difference signal map via a — Ds/AD^ 1 . Considering all the angular scales in 
a CMB experiment like Planck 2 they find a — 2 after t — 2 ■ 10 yr. Performing a harmonic transform of the map they 
then isolate the £ = 1 mode and find for the same time span a — 40 with AD$ = 3AD^ / (Air). 

Using just the information on the variation of the power spectrum SCi, where signal-to-noise ratio is a — 
\8C\\/2ACi , gives a lower statistical significance (a = 4) due to the removal of the time covariance information. As 
already mentioned previously, the variance of the C/s 

AC? = y^ T (2<C| st > N +N e ) (105) 

does not contain the cosmic variance term, but includes the uncertainty on the estimated power spectrum in a certain 
noise realisation ensemble. In Fig. 22 a ACDM temperature anisotropy power spectrum is compared to its variance 
(neglecting the intrinsic cosmic variance) and the simple noise y/2/ f s k y /(2£ + TjNf for a Planck-like experiment. The 
difference between the two noises arises because the power spectrum is estimated from a map: after a certain number 



2 http://www.esa.int/planck 
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Figure 22. CMB temperature anisotropy power spectrum (solid curve) confronted with the error ACf (dashed curve) and the 
pure noise s/2/f sky /{2l+l)Nt (dot-dashed curve). 



of noise realisation the expected value of the spectrum is < Cf st > with a variance that takes into account both the 
noise and the scatter in the estimate [71]. 

Regarding the dipole, the left panel of Fig. 23 shows the signal-to-noise ratio relative to it for a Planck-like mission 
and a 100-fold enhanced mission. Using the difference map and after removing the small scale contamination it should 
be possible to detect a time variation of the dipole in lOyr. 

Following the formalism in [81] the authors in [13] calculated the difference in high precision maps taken a century 
apart and compared it to the experimental uncertainty, that is the noise with 

Ni = ^J% hm exp (£ 2 ^, im /8m2), (106) 

where 0f w h m is the FWHM of the gaussian beam profile. The predicted experiment has a sensitivity of s = 40mK 
s 1 / 2 (where a p i X — s/y/t p i X ). an angular resolution of 0.86', a 4 year time mission covering 75% of the sky. The right 
panel of Fig. 23 shows the power spectrum of the difference of two maps with a time interval of 100 yr, together with 
the errors after a binning in £ and the noise curve of such an experiment . 

Among all the contamination that might smudge the detection, the uncertainty on the calibration seems to be the 
less alarming due to the fact that we are dealing with a difference map. The less promising systematic is probably 
caused by the time-varying foreground emission. The observer foreground might change with the same time scale as 
the dipole. Hopefully, multi- wavelength observations will help constraining them. However, it would take thousands 
of years to detect a change in the spectrum at higher multipoles. 



VIII. CONCLUSIONS 



Time domain information provides critical insight in many branches of astrophysics, even cosmology (distant su- 
pernovae and active galactic nuclei being obvious examples). However, the universe as a whole is evolving slowly 
and so one expects this evolution to be difficult to observe. Real-time cosmology tests directly cosmic kinematics by 
observing changes in source positions and velocities. 

Measurements of the redshift drift (or velocity shift) with future extremely large telescopes and high-resolution 
spectrographs could provide interesting information on the source of cosmic acceleration, which would complement 
other, more traditional cosmological tools. The observation of redshift drift alone can be affected by strong param- 
eter degeneracies, limiting its ability to constrain cosmological models. Uncertainties on parameter reconstruction 
(particularly for non-standard dark energy models with many parameters) can be rather large unless strong external 
priors are assumed. When combined with external inputs, however, the time evolution of redshift could discriminate 
among otherwise indistinguishable models. 
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Figure 23. Left panel: Signal-to-noise ratio as a function of time for observing a change in the CMB dipole. Thin lines are for 
a Planck-like instrument and thick lines are for an instrument 100 times as sensitive. Solid curves are drawn from differences 
in maps, while dashed curved after removing small scale noise and dotted curves contain no time-time correlation information. 
From Ref. [71]. Right panel: The bottom curve is the power spectrum of the difference between two maps taken a century apart. 
The dotted line corresponds to the noise (Eq. (106)) and the hashed boxes are the errors after binning in I. For comparison 
the authors superposed on the same plot the EE power spectrum for r = 0.1, the BB power spectrum for two different values 
of the tensor-to-scalar ratio. From Ref. [13] 



Despite its inherent difficulties, the method has many interesting advantages. One is that it is a direct probe of 
the dynamics of the expansion, while other tools (e.g. those based on the angular diameter distance) are essentially 
geometrical in nature. This could shed some light on the physical mechanism driving the acceleration. For example, 
even if the accuracy of future measurements turns out to be insufficient to discriminate among specific models, this 
test would be still valuable as a tool to support the accelerated expansion in an independent way, or to check the 
dynamical behaviour of the expansion expected in general relativity compared to alternative scenarios. In LTB void 
models for example, its constraining power is particularly evident since the LTB expansion is always decelerated and 
the effect turns out to be mostly sensitive to the scale of the void, but not to other particular void properties like 
steepness of the transition. 

Furthermore, despite being observationally challenging, the method is conceptually extremely simple. For example, 
it does not rely on the calibration of standard candles (as it is the case of type la SNIa) or on a standard ruler 
which originates from the growth of perturbations (such as the acoustic scale for the CMB or for transverse BAO) 
and on effects that depend on the clustering of matter (except on scales where peculiar accelerations start to play a 
significant role). Finally, it is at least conceivable that suitable sources at lower redshifts could be used to monitor 
the redshift drift in the future. This would be extremely valuable, since some non-standard models have a stronger 
parameter dependence at low and intermediate redshifts, that could be exploited as a discriminating tool. Exploring 
the feasibility of such proposals will certainly be an interesting topic for further studies from observers. 

Peculiar and proper accelerations would also be an interesting target for such future measurements, since they 
could give an independent measurement of the mass profile of clusters of galaxies, and possibly even of individual 
galaxies. Measurements of the velocity shift in real time for different mass scales could test the relation between the 
concentration parameter of clusters and the virial mass. 

The expected signal for some objects (like Coma cluster) might be comparable to the cosmological change in redshift 
expected for distant quasars due to the expansion dynamics. It then seems plausible that it could be observable as 
well, although the objects used as test particles in the two cases are different in nature. This would also be useful to 
the extent of assessing the effect of peculiar acceleration as a source of errors in the measurement of the cosmological 
velocity shift. 

Experimental and observational details of future instruments are still too blurry to make any clear statement on 
the actual feasibility of using peculiar acceleration as an astrophysical and cosmological tool (even less for proper 
acceleration). Many issues are still to be discussed before forecasting the results of actual observations, from the 
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stability of the galactic spectra to the effect of local deviations from spherical symmetry. However, it is exciting to 
entertain the possibility that peculiar and also proper accelerations might actually be observable in the not too distant 
future, and hope that more detailed investigations of their applicability in astrophysics might take place in the next 
few years. 

In particular, the possibility of reconstructing the gravitational potential of a galaxy by means of the redshift drift 
signal is also appealing. Assuming that either the cosmological redshift of the galaxy is almost vanishing or it can be 
averaged out, the velocity shift produced by peculiar acceleration of test particles orbiting outside the galactic disc 
might be a probe of a modified gravity configuration, in which the CDM halo is absent and the Poisson equation 
is modified following MOND: although the two scenarios are undistinguishable using rotation velocity curve, the 
question of whether peculiar acceleration might be employed to discriminate between them seems quite interesting. 

As expected, the closer is the test particle to the galactic center the higher is the signal, reaching values compara- 
ble to the cosmological ones predicted by simulations and possibly detectable from high-resolution and ultra-stable 
spectrograph coupled to new generation telescope. 

In our Galaxy, the 150 known globular clusters orbiting outside the disc could play the role of test particles. The 
advantage of adopting globular cluster as test particles is related to the fact that they not only are bright, but more 
importantly that they are made of thousands of stars and that their position is fairly well known. As in this case 
we clearly are internal observers, the projected acceleration of the Sun on the disc, which has its maximum value on 
the order of a few cm/s, is assumed to be subtracted, and the cluster-centered acceleration of single stars in globular 
clusters is assumed to be averaged out. The effect of being off-centered observers combined with the non-radial 
direction of the acceleration in Kuzmin disc makes the pattern of the contours non trivial. However, while the signal 
reaches the maximum value close to the Galactic centre for both scenarios, the difference signal is stronger at high 
Galactic longitude where the spherical CDM halo is more influent. 

Although the feasibility of such observations is still under debate, this new technique would open a new window to 
the ability of testing the gravitational potential and of constraining the parameters involved in its expression, such as 
the mass of the galaxy. 

Any anisotropy will leave an imprint on the angular distribution of objects that are able to trace cosmic expansion. 
Planned space-based astrometric missions aim at accuracies of the order of a few micro-arcseconds. In the framework of 
real-time cosmology, the cosmic parallax of distant sources in an anisotropic expansion might be observable employing 
the same missions. A positive detection of large-scale cosmic parallax would disprove therefore one of the basic tenets 
of modern cosmology, isotropy. This test may probe a different range of scales depending on the quasar redshift 
distribution. Contrary to CMB limits, the cosmic parallax method cannot be completely undermined by the observer's 
peculiar motion, and is limited only by source statistics instead of cosmic variance. The major source of systematics 
would be the subtraction of the aberration change parallax. 

If an anisotropy is present before the last scattering surface, CMB maps will also be affected. The temperature field 
will carry extra anisotropics mainly caused by the angular dependence of the redshift at decoupling. By resolving 
geodesic equations and expanding temperature anisotropics in spherical harmonics, it is straightforward to relate the 
low multipole components to the eccentricities of the model. In Bianchi I models the first notable multipoles related 
to the CMB are the monopole and the quadrupole. The observed value of the latter puts constraints on the shear 
at last scattering of order 10 -5 , taking into account the cosmic variance. In this still exploratory stage of analysis 
it is worth stressing that CMB puts indeed a very powerful constraint on the shear at the time of decoupling, but 
with almost no direct impact on late time expansion history. Complementary to that, cosmic parallax is a direct and 
potentially powerful test of anisotropy at small redshifts. 

The anisotropic stress of dark energy is expected to have a leading role in the generation of anisotropy at late times. 
It can be parameterized by skewness parameters in the stress-energy tensor formulations, which may be constant or 
time dependent functions. The only way to test these models is to use either the angular dependence of the magnitude 
or the angular distribution of objects in the sky at recent time, i.e. either distant source angular distribution or the 
real-time cosmic parallax, the two relying on different techniques and having independent systematics that complement 
each other. 

Gaia will be able to constrain the skewness parameters up to 10~ 3 -j- 10~ 4 at 2a, comparable to CMB tests at 
decoupling time (a Gaia+ experiment would improve them by one order of magnitude). 

Unlike the LTB models with off-centre observers, the cosmic parallax signal in Bianchi I models is a combination 
of two quadrupole functions of the two angular coordinates. Since the most important systematic noises, caused 
by peculiar velocities and aberration changes, have a dipolar functional form, Bianchi I models seem to be ideally 
testable, though even in LTB models specific observational strategy aiming at distinguishing the signal from the noise 
are possible. 

CMB and cosmic parallax detect anisotropy at two different times and, from an observational point of view, are 
completely independent on each other: by combining them together one will have the opportunity to reconstruct the 
evolution of the anisotropy and test with high accuracy the Copernican Principle. 
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The two effects, the redshift drift together with cosmic parallax, can be detected with the E-ELT and with Gaia 
or an enhanced version of Gaia. These tests can be employed to distinguish a LTB void from an accelerating FRW 
universe, possibly eliminating an exotic alternative explanation to dark energy. 

A 4cr separation can be achieved with E-ELT in less than 10 years, much before the same experiment will be able 
to distinguish between competing models of dark energy. A Gaia-like mission, on the other hand, can only achieve a 
reasonable detection of a void-induced cosmic parallax in the course of 30 years. 

Nevertheless, cosmic parallax remains an important tool and in fact one of the most promising way to probe 
general late-time cosmological anisotropy. In particular, even if it only lasts 6 years, Gaia should constrain late-time 
anisotropics similarly to current supernovae catalogs, but in an independent way. Also, in a FRW model it can be 
used to measure our own peculiar velocity with respect to the quasar reference frame and consequently to the CMB, 
therefore providing a new and promising way to break the degeneracy between the intrinsic CMB dipole and our own 
peculiar velocity. Combined, they will form an important direct test of the FRW metric. 

As any proper motion signal increases linearly with time, any future mission with a global astrometric accuracy as 
good as Gaia can be used to detect the cosmic parallax signal. 

It's really compelling that two great tools like Gaia and E-ELT are being planned now when we begin to realize 
the importance of ultra-precise astrometric and spectroscopic measurements for cosmology. 

The question of how long one would have to wait to see new information from the CMB has also been addressed. 
We expect to observe a change on time scales much shorter than cosmological. By considering optimal estimates 
for differences in observed skies, using currently available detector technology like Planck, it turns out that the 
dipole might be observed to change due to our Galactic motion in about a decade, while unfortunately it would take 
thousands of years to detect even a statistical change in the higher order multipoles. 
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Appendix A: Geodesies in LTB Models 

Due to the axial symmetry and the fact that photons follow a path which preserves the 4- velocity identity u a u a = 0, 
the four second-order geodesic equations for (t, r, 9, if) can be written as five first-order ones. We will choose as variables 
the center-based coordinates t, r, 9, p = dr/dX and the redshift z, where A is the affine parameter of the geodesies. 
We shall refer also to the conserved angular momentum 



J = R 2 — = const = J Q . (Al) 
dA 

For a particular source, the angle £ is the coordinate equivalent to 6 for the observer, and in particular £ is the 
coordinate £ of a photon that arrives at the observer at the time of observation to. Obviously this coincides with the 
measured position in the sky of such a source at to. 

As per Figure 6 and Section II B 1 we will refer to (t, r, 9, (f) as the comoving coordinates with origin on the center 
of a spherically symmetric model. Peculiar velocities aside, the symmetry of such a model forces objects to expand 
radially outwards, keeping r. 9 and (f> constant. In terms of these variables, and defining A such that u(A) < 0, the 
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autonomous system governing the geodesies is written as 



dt 
dA 
dr 
dA 
d0_ 
dA 
dz 
dA 

dp 
dA 




P, 
J 



(A2) 



(1 + *) 



R' R 9 R T 9 
— P + R* J 




2R'p 



2 + 2/3 



J?' 



The angle £ along a geodesic is given by [38] : 



C0S C = — 



R'(t,r)p 



(A3) 



from which we obtain, exploiting the remaining freedom in the definition of A, the relations [38] Jo = J = 
R(to, ro) sin(£o) and po = — yl + /3(?"o) cos(£o) / R'(to, r o)- Therefore, our autonomous system is completely de- 
fined by the initial conditions to, ro, Oq — 0, zq — and £o- The first two define the instant of measurement and the 
offset between observer and center, while £ stands for the direction of incidence of the photons. 

Following [6] , an algorithm for predicting the variation of an arbitrary angular separation and redshift with time 
can be written as follows: 

1. Denote with (z a i,£ a i) the observed coordinate of a source at a given time to and observer position ro; 

2. Solve numerically the autonomous system with initial conditions (to,ro,0o = 0, zo = 0,£o = £ai) and find out 
the values of A* such that z(\*) = z a \ ; 

3. Take note of the values r a i(A*) and Q i(A*) (since the sources are assumed comoving with no peculiar velocities, 
these values are constant in time); 

4. Define as the parameter value for which r a2 (Xl L ) — r al (A*), where r a2 is the geodesic solution for a photon 
arriving a time At later with an incident angle £ a 2, and vary £ a 2 until 6* a 2 ( Ajj) = #ai(A*) ; 

5. Compute the difference Az = z a i — z a \ (to obtain the redshift drift); 

(£02 - 62) - (ia\ - ib\) (to 



6. Repeat the above steps for source b, and compute the difference A t j 
obtain the cosmic parallax). 



72 - 7i 



The above algorithm provide two (in principle) coupled observables: the Cosmic Parallax (CP) (see Sec. Ill), and 
the redshift drift (see Sec. II B 1). 

A remark on the above procedure is in order. Due to the intrinsically smallness of both the cosmic parallax and 
the redshift drift (in the course of a decade), a carefully constructed numerical code is needed to correctly compute 
either. To give an idea of the amount of precision required, consider the following: if one naively calculates A t 7 for 
a At of 10 years, one needs to evaluate £ a i and £ a 2 with at least 13 digits of precision (as the CP is of the order of 
0.2 ^as - 10- 12 rad) [25]. 



Appendix B: Error Bars Estimates with SDSS Quasars 

The signal-to-noise ratio per pixel (of size 0.0125 A) in (103) was estimated in [9] to be 



D 

42 m 



tint £ 

lOh 0/25 



(Bl) 



where Zx and mx are the source zero point and apparent magnitude in the "X" band and D, t int and e are the 
telescope diameter, total integration time and total efficiency respectively. We assumed a central obscuration of the 
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telescope's primary collecting area of 10% [9]. Note that D = 42 m corresponds to the reference design for the 
E-ELT [82]. 

The reason magnitudes are quoted in terms of an arbitrary "X" band is because one should use the magnitude of the 
bluest filter that still lies entirely redwards of the quasar's Lyman-a emission line [9]. This means that for zqso < 2.2 
one should use the magnitude in the g-band; for 2.2 < zqso < 3.47 the one in the r-band; for 3.47 < zqso < 4.61 
the i-band; for zqso > 4.61 the z-band. A good estimate for mx can be achieved with the SDSS DR7, selecting the 
brightest quasars in each redshift bin using the appropriate band for such bin. 

Following [22], Ref. [25] selected 40 quasars in 5 redshift bins, centered at z = {2, 2.75, 3.5, 4.25, 5}, all of the same 
redshift width of 0.75. The corresponding bands are, in order, {g, r, r, i, z} (where the i-band could equally be chosen 
for the middle bin). Doing so, one gets for the average (amongst the 8 brightest quasars) apparent magnitude mx 
for each bin the following: mx = {15.45, 16.54, 16.40, 17.51, 18.33}. Finally, the zero point magnitude ratio in each 
bin was estimated in [25] to be: Zx/Z r = {1.01, 1.00, 1.00, 0.98, 0.93}. The accuracy of this last estimate is however 
quite unimportant in the results shown. 
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